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Abstract 

Presented  in  this  work  are  the  results  of  an  investigation  of  alternative  means  for 
powering  spacecraft  and  launch  vehicles  with  energy  sources  other  than  chemical 
combustion.  Nuclear  Thermal  Propulsion  (NTR)  and  the  energy  release  of  a  nuclear  spin 
isomer  present  potential  for  increased  rocket  performance  with  a  compact,  high-energy 
fuel  sources  replacing  the  combustion  engines  of  the  Delta  IV-H  1st  and  2nd  stage 
vehicles.  NTR  was  represented  by  the  Nuclear  Engine  for  Rocket  Vehicle  Application, 
CERMET,  and  the  Particle  Bed  Reactor  (PBR)  fission  designs,  while  the  isomer 
hafnium-  178-m2  was  investigated  in  a  PBR  configuration.  Energy  storage  levels  of  1.3 
GJ/g  are  possible  with  this  material,  though  the  successful  triggering  and  maintenance  of 
a  chain  reaction  in  this  material  are  still  debated  topics  within  the  scientific  community. 
The  best  application  for  either  technology  is  as  an  upper  stage  vehicle  with  the  shielding 
requirements  reduced  to  that  of  just  a  shadow  shield  between  the  core  and  the  spacecrafts 
upper  structure.  The  fission  designs  are  capable  of  specific  impulse  values  between  800 
and  1,000  s  leading  to  mass  savings  in  the  range  of  7,000  to  nearly  10,000  kg  once  the 
engine  masses  and  shielding  have  been  included.  An  isomer  core  in  the  configuration  of 
a  19-element  PBR  may  be  able  to  achieve  a  specific  impulse  on  the  order  of  880  s  with 
the  isomer  in  metallic  fonn,  and  specific  impulse  values  as  high  as  1 ,090  s  if  the  isomer  is 
in  the  fonn  of  hafnium  carbide.  This  translates  to  somewhere  between  a  5,000  and  9,000 
kg  depending  on  the  material  makeup  of  the  core  and  heat  efficiency.  Payload  mass 
increases  by  a  factor  of  two  or  greater  velocity  change  capability  are  the  payoffs  of  these 
systems. 
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ISOMER  ENERGY  SOURCE  FOR  SPACE  PROPULSION  SYSTEMS 


I.  Introduction 


1.1.  Motivation 

The  current  capabilities  of  our  space  propulsions  systems  are  limiting  human 
operation  within  and  about  the  space  environment.  Chemical  rocket  systems,  while 
necessary  for  launch  operations  due  to  their  high  thrust  characteristics,  do  not  achieve  the 
highest  levels  of  fuel  economy.  This  economy  is  measured  by  a  vehicle’s  specific 
impulse  (Isp)  value  which  represents  the  time  a  rocket  can  produce  one  pound  of  thrust 
using  one  pound  of  fuel.  Further  complexities  within  a  rocket  design  are  introduced  with 
the  choice  of  relying  on  chemical  propulsion;  since  both  a  fuel  and  oxidizer  are  required, 
must  be  stored  separately,  and  be  fed  to  the  combustion  chamber.  Near  earth  propulsion 
systems  relying  on  solar  power,  while  efficient  in  terms  of  fuel  economy  (solar  thennal) 
or  requiring  no  fuel  at  all  (solar- lightsail),  lack  the  thrust  required  for  rapid  acceleration 
and  short  time-scale  missions.  The  inability  of  chemical  and  solar  powered  spacecraft  to 
both  “rapidly  and  efficiently”  operate  in  the  space  enviromnent  is  delaying  mankind’s 
progress  in  its  quest  to  explore,  utilize  space  assets  to  better  life  on  Earth,  and  expand  our 
civilization  in  such  a  manner  that  it  allows  us  to  grow  as  a  nation  and  world  community. 
The  new  challenges  presented  by  life  in  space  will  lead  to  new  ideas,  new  technologies, 
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and  ultimately  to  an  evolution  of  our  society.  While  there  will  always  be  those  among  us 
content  with  living  in  a  society  that  does  not  extend  its  reach  into  space,  many  people 
now  realize  the  importance  of  space  and  the  benefits  that  conquering  this  domain  will 
have  on  our  way  of  life  (1).  Settling  this  new  frontier,  or  establishing  a  pennanent 
presence  in  space  is  important,  and  it  is  rapidly  being  realized  that  the  use  of  energy 
sources  other  than  chemical  or  solar  are  needed  to  make  this  vision  a  reality. 

While  it  is  possible  to  operate  near-earth  with  chemical  or  solar  powered  systems, 
the  key  to  human  advancement  is  efficient  and  rapid  transport  reducing  the  timescales  of 
operations  in  the  space  environment.  What  is  required  is  an  alternative  fuel  source  high 
in  energy  density  and  capable  of  releasing  this  energy  on  demand.  For  the  past  sixty 
years  scientists  and  engineers  have  contemplated,  researched,  engineered,  designed, 
tested,  and  even  flown  a  few  systems  powered  by  nuclear  fission  or  radioactive  decay 
(2:10-15;  3:452-459;  4:12-19;  5;  6).  These  vehicles  ranged  from  launch  systems  to  deep 
space  probes  and  spacecraft. 

Deep  space  missions  are  most  often  discussed  as  the  target  for  space  nuclear 
reactors;  this  is  more  an  area  of  interest  for  the  National  Air  and  Space  Administration 
(NASA)  (7)  than  for  the  Air  Force,  and  not  the  only  need  for  such  propulsion  sources. 
Near  tenn  there  is  a  need  for  increasing  payload  amounts  delivered  to  space  and  for 
missions  reaching  out  as  far  as  the  moon  (8),  therefore;  an  investigation  of  nuclear 
electric  propulsion  was  not  included  in  this  study. 

While  concepts  like  controlled  nuclear  fusion  and  antimatter  systems  are  still  in 
the  far-off  future,  fission  reactors  (well  understood  and  operated  for  terrestrial  energy 
generation  and  naval  vessel  propulsion)  would  appear  to  be  the  near-tenn  solution  to  the 
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problem.  However,  many  hazards  are  associated  with  the  production,  storage,  handling, 
and  use  of  fission  reactors,  which  is  why  today  we  are  still  relying  on  chemical  systems. 
If  the  hazards  of  fission  reactors  cannot  be  overcome  or  accepted,  another  source  of 
energy  is  required  to  satisfy  this  need.  It  is  possible  to  achieve  metastable  excited  energy 
states  in  some  nuclei  which  can  serve  as  a  means  of  energy  storage.  Nuclear  spin 
isomers,  in  particular,  the  isomer  hafnium- 178- m2  (178Hfn2)  stores  approximately  2.446 
MeV  per  atom  or  1.3  GJ/g  with  a  31 -year  half-life  (9).  Research  and  experiments 
conducted  over  the  past  5  years  have  indicated  that  it  may  be  possible  to  trigger  this 
isomer  to  release  its  stored  energy  on  demand  (10:4).  This  is  indeed  remarkable  since 
gram  quantities  can  store  energies  equivalent  to  metric  ton  quantities  of  chemical  fuels  or 
explosives  (10:4).  In  addition,  the  spectrum  of  ionizing  radiation  released  by  this  decay 
is  different  from  that  seen  in  a  conventional  fission  reaction.  Nuclear  fission  results  in 
the  production  of  fission  fragments,  neutrons,  alpha  and  beta  particles,  and  gamma-rays. 
Gamma-rays,  in  particular,  can  range  from  0.2  to  7.6  MeV  in  energy,  from  the  fission  of 
U235  (1 1:7-72).  With  the  decay  of  178Hfn2  the  energy  released  is  entirely  in  the  form  of 
gamma-rays  ranging  from  12.7  to  547  keV.  If  this  released  energy  could  be  used  to  heat 
a  propellant  (directly  or  indirectly)  to  significant  temperatures  then  there  is  potential  for  a 
new  source  of  energy  to  transport  us  rapidly  through  the  space  environment. 

1.2.  Research  Objective 

The  goal  of  this  research  was  to  detennine  how  triggered  isomer  decay  could  best 
be  utilized  as  a  source  of  heat  energy  in  a  nuclear  thermal  rocket  (NTR)  configuration 
This  includes  study  of  the  current  state  of  triggered  isomer  research  an  investigation  of 
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the  interaction  of  isomer  decay  products  with  materials,  research  into  the  best  suited 
nuclear  reactor  design,  and  comparisons  against  chemical  and  nuclear  propulsion  systems 
as  they  are  used  in  fielded  systems  today  or  proposed  to  be  used  in  the  systems  of 
tomorrow. 

The  successful  triggering  of  the  isomer  1 7SHfl2  is  a  highly- debated  topic.  To  date 
only  a  few  experiments  have  demonstrated  successful  triggering  of  this  material,  and 
these  have  not  been  universally  accepted  by  the  scientific  community(12;  13;  14;  15;  16; 
17).  The  belief  that  triggering  is  possible  is  necessary  for  this  research,  but  a  summary  of 
the  research  conducted  by  scientists  on  both  side  of  the  debate  is  included  in  Sec.  2.1  of 
this  document  to  accurately  represent  the  maturity  of  triggered  isomer  technology. 

All  debates  about  the  triggering  of  178Hfn2  aside,  the  primary  concern  becomes 
absorbing  the  electromagnetic  radiation  released  by  the  radionuclide  decay.  The 
attenuation,  and  energy  deposited  by  this  gamma  radiation  becomes  important  when  one 
tries  to  heat  a  propellant  or  shield  a  spacecraft  and  its  surroundings  from  radiation  given 
off  by  the  core. 

The  most  heavily  researched  thermal  reactors  designed  for  space  use  are  the 
Nuclear  Engine  for  Rocket  Vehicle  Application  (NERVA),  the  Particle  Bed  Reactor 
(PBR),  and  the  CERMET  reactor  (3:452-459;  4:12-19;  5;  6;  18;  19;  20;  21).  While  none 
of  these  reactors  ever  made  it  to  production,  enough  testing  was  conducted  to  prove  the 
concepts,  especially  in  the  case  of  the  NERVA  rocket  which  made  it  all  the  way  to  full- 
scale  testing.  Table  1  lists  the  capabilities  of  these  reactors  as  they  were  designed.  Each 
reactor  design  is  unique  in  its  operating  principal  but  this  study  will  aim  at  determining 
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the  traits  of  each  reactor  design  that  would  be  most  beneficial  to  a  triggered  isomer 
design. 


Table  1.  Comparison  of  Possible  Near-Term  Concepts  for  Reactors  (3:457) 


Nuclear  Engine 

NERVA 

CERMET 

PARTICLE- 

BED 

Power  (MW) 

1,570 

2,000 

1,945 

Thrust  (N) 

334,061 

445,267 

333,617 

Propellant 

h2 

h2 

h2 

Fuel  Element 

Solid  rod 

Solid  rod 

Porous  particle 
bed 

Maximum  Propellant 
Temperature  (K) 

2,361 

2,507 

3,200 

Isp  (s) 

825 

930 

971 

m  (kg/s) 

41.27 

48.81 

35.02 

Chamber  Pressure 
(MPa) 

3.102 

4.136 

6.176 

Nozzle  Expansion 
Ratio 

100 

120 

125 

Engine  Mass  (kg) 

10,138 

9,091 

1,705 

Total  Shield  Mass 
(kg) 

1,590 

1,590 

1,590 

Engine  F/W  ( no 
shield) 

3.4 

5.0 

20.0 

A  similar  core  design  replacing  the  fission  fuel  with  the  isomer  1 7sHlm2  is  the  starting 
point  for  this  work.  This  will  have  to  be  carried  out  in  such  a  manner  that  a  chain 
reaction  is  ignitable  and  maintainable,  though  it  is  realized  at  this  time  that  scientists  and 
engineers  may  be  a  long  way  from  achieving  this  experimentally.  In  essence,  the  best 
reactor  design  to  support  triggered  isomer  decay  will  be  identified  through  an  analysis  of 
materials  and  geometry.  This  will  be  followed  by  a  comparison  against  the  baseline 
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fission  reactor  and  then  comparisons  against  chemical  systems  used  in  the  selected 
mission  to  see  if  there  is  a  potential  for  perfonnance  improvements. 

The  vehicle  that  was  examined  in  this  study  was  the  Delta  IV-Heavy  launch 
vehicle,  both  its  1st  and  2nd  stages.  The  Delta  IV-H  rocket  was  chosen  because  it  is  one  of 
the  most  powerful  rockets  that  will  be  used  to  boost  payloads  into  orbit  in  the  next  few 
years.  It  also  already  incorporates  cryogenically-  fueled  stages  that  hopefully  will 
eliminate  the  need  for  drastic  alterations  to  the  propellant  storage  and  handling  system 
with  an  alternative  design.  Table  2  lists  the  perfonnance  characteristics  of  each  of  these 
stages. 


Table  2.  Vehicle  Design  Parameters  (22) 


Vehicle  Name 

1st  Stage  Delta  IV-H 

2na  Stage  Delta  IV-H 

Inert  Mass  (kg) 

80,280 

3,490 

Propellant  Mass  (kg) 

598,800 

27,200 

Vacuum  Thrust  (N) 

9,945,000 
{3  engines} 

110,000 
{1  engine} 

Chamber  Pressure 
(MPa) 

9.72 

3.21 

Vacuum  Isp  (s) 

420.0 

462.4 

Nozzle  Expansion  Ratio 

21.5:1 

285:1 

Propellant  Mass 
Fraction 

0.88 

0.89 

Vehicle  F/W( with 

10,843  kg  payload) 

1.4 

0.27 

Av  (km/s) 

7.84* 

4.62* 

The  thrust  values  shown  are  average  vacuum  thrust  levels  and  the  number  in  brackets 
indicate  the  combined  number  of  engines  producing  this  thrust.  The  velocity  values 
(asterisked)  were  calculated  based  on  specific  mission  goals  and  these  calculations  are 
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displayed  in  Appendix  F.  In  the  case  of  the  Delta  IV-H  1st  stage  the  velocity  change 
(Av  )  calculated  was  that  velocity  needed  to  enter  a  circular  orbit  at  an  altitude  of  1 10  km 
above  the  Earth.  The  2nd  stage  velocity  change  was  determined  from  knowledge  of  the 
2nd  stage  specific  impulse,  burn  time,  vehicle  mass,  and  thrust  value  assuming  constant 
thrust.  This  allowed  for  detennining  the  mass  flow  rate,  propellant  exhaust  velocity,  and 
finally  its  Av  through  use  of  the  ideal  rocket  equation  shown  in  Eq.  [23].  Once  a  reactor 
design  is  identified  that  has  the  highest  potential  for  successfully  operating  on  the 
principals  of  triggered  isomer  decay,  the  performance  will  be  compared  against  what  is 
predicted  for  the  current  design  which  is  combustion  driven 

1.3.  Past  Work. 

Prior  to  this  research  effort  two  other  studies,  both  by  master’s  students  at  the  Air 
Force  Institute  of  Technology,  focused  on  propulsion  applications  for  using  triggered 
isomer  energy.  The  first  study  was  conducted  by  Captain  Carl  Hartsfield,  and  was 
entitled  “Analysis  of  the  Application  of  a  Triggered  Isomer  Heat  Exchanger  as  a 
Replacement  for  the  Combustion  Chamber  in  an  Off-the-Shelf  Turbojet”  (23),  completed 
in  March  2001.  Hartsfield  modeled  variations  of  flat  plate  configuration,  solid-state  heat 
exchangers  with  the  commercial  software  package,  ANSYS®  5.6. 1  (24).  Three  basic  heat 
exchanger  geometries  were  studied,  and  Hartsfield  concluded  that  all  could  produce 
sufficient  heat  transfer  to  replace  the  combustion  chamber  in  the  J-57  turbojet  engine 
used  in  the  Boeing  707,  KC- 135,  and  B-52  (23  :XV).  In  this  design  the  heat  exchanger 
was  assumed  to  be  manufactured  from  the  isomer  material  and  heater  exit  temperatures 
resulted  in  the  range  of 986  to  1,150  K  with  sea-level  thrust  values  from  37,000  to  47,000 
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N.  Though  constant  heat  generation  is  likely  the  most  accurate  model  of  the  isomer 
decay  process,  it  was  deemed  too  computationally  complex  and  a  constant  surface 
temperature  of  2,400  K  was  applied  to  all  the  surfaces  of  the  heat  exchanger.  This 
decision  between  modeling  with  constant  surface  temperature  and  constant  surface  heat 
flux  was  made  on  the  basis  of  the  temperature  gradient  for  the  constant  surface 
temperature  condition  most  nearly  matching  that  of  the  constant  heat  generation 
condition.  All  of  the  radiation  escaping  the  engine  was  assumed  to  be  in  the  form  of  600 
keV  gamma-rays,  which  is  slightly  higher  than  the  most  energetic  photon  emitted  during 
the  natural  decay  of  178Hfn2  which  is  574  keV  (25).  Only  5%  of  the  heat  generated  in  the 
decay  process  was  assumed  to  escape  the  heat  exchanger  by  radiation  to  the 
surroundings.  Stagnation  temperatures,  pressure  losses,  radiation  shielding,  and  the 
physical  dimension  of  a  triggered  isomer  heat  exchanger  were  all  examined  in  this  study 
with  the  result  being  a  predicted  heat  exchanger  volume  of  0.042  m3  and  weight  of 
approximately  420  kg. 

Captain  Chris  Hamilton  followed  up  this  work  with  a  “Design  Study  of  Triggered 
Isomer  Heat  Exchanger- Combustion  Hybrid  Jet  Engine  for  High  Altitude  Flight”  (2),  in 
March  2002.  A  Global  Hawk  Unmanned  Aerial  Vehicle  (UAV)  relying  on  a 
conventional  combustion  engine  for  takeoff  and  near-earth  operation,  and  a  triggered 
isomer  heat  source  for  high  altitude  flight  (>20,000  ft)  was  the  subject  of  this  study.  It 
was  concluded  in  this  work  that  a  single  hybrid  engine  with  a  switchover  from 
conventional  combustion  to  Triggered  Isomer  Heat  Exchanger  (TIHE)  operation  could 
extend  the  endurance  of  a  High  Altitude  Long  Endurance  Aircraft  (HALE)  weeks  beyond 
the  current  mission  limitation  which  was  on  the  order  of  days  (2).  Up  to  a  20%  drop  in 
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vehicle  weight  was  also  possible  primarily  due  to  lower  fuel  requirements  (2).  Weight 
estimates  including  8.9  kN  for  the  triggering  photon  source  and  engine  modifications, 

1 1.6  kN  for  semispherical  shielding,  and  6.67  kN  for  the  TIHE,  were  also  made  (2). 
Hamilton  used  Aircraft  Engine  Design  System  Analysis  Software  (AEDsys),  version  2.13 
(26)  and  On-Design  Analysis  of  Gas  Turbine  Engines  (ONX),  version  4.021  (27)  for 
basic  engine  design  and  performance  analysis.  The  same  assumptions  from  Hartsfield’s 
study  concerning  the  heat  exchanger  design  were  carried  over  into  this  study  with  the 
addition  of  the  requirement  that  the  heat  exchanger  produce  equal  heating  rates  to  the 
chemical  combustors  it  was  replacing.  The  development  of  a  specific  reactor  design  and 
its  method  of  manufacture  were  left  as  areas  for  further  investigation  by  both  studies.  No 
other  studies  of  using  a  triggered  isomer  energy  source  for  propulsion  were  found  in  the 
literature  search  performed. 
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2.  Background/History 


2.1.  Triggered  Isomer  Development 

The  4- quasiparticle  isomer  of  hafnium  with  excitation  energy  equal  to  2.446  MeV 
and  a  31-yr  half- life  (9)  represents  a  compact,  long  duration  means  of  energy  storage  that 
recent  experiments  (9;  10;  28;  29;  30:2-3)  have  shown  can  be  manipulated  to  release  its 
energy  on  demand.  Some  shell- model  calculations  have  also  suggested  the  existence  of 
more  extreme  isomers  such  as  the  neutron- rich  hafnium  isotopes,  and  in  the  isotopes  of 
lutetium  (Z  =  71)  and  tantalum  (Z  =  73)  (3 1 :84-85),  but  the  amount  of  research  put  into 
the  study  of  the  hafnium  isomer  to  date  makes  it  the  most  likely  candidate  for  actual 
applications.  The  tenn  nuclear  isomer  is  used  to  describe  long-lived,  high-energy  states 
of  excitation  in  nuclei  as  opposed  to  the  term  chemical  isomer  which  describes  a 
variation  in  the  bond  arrangement  of  a  molecule.  The  hafnium  nuclear  isomer  can  be 
fonned  by  bombarding  tantalum  with  protons  in  a  process  that  requires  a  nuclear  reactor 
or  particle  accelerator.  Quantities  of  this  isomer  are  also  produced  in  Dubna,  Russia  by 
bombarding  176Yb  with  alpha-particles.  Bombardment  of  these  materials  results  in  decay 
to  the  l78Hf12  isomer,  in  the  instance  where  energy  is  transferred  to  the  excited  states  of 
product  nuclei  instead  of  manifesting  as  kinetic  energy  of  escaping  products.  At  present 
only  small  quantities  of  the  isomer  can  be  produced  and  are  available  for 
experimentation.  SRS  Technologies  in  Huntsville,  Alabama,  is  under  contract  with  the 
Air  Force  Research  Laboratory  for  a  supply  of  the  isomer,  but  at  present  only  a  ten- 
thousandth  of  a  gram  quantity  can  be  generated  at  a  given  time  and  the  cost  of  production 
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is  likely  to  be  expensive  (32).  A  near-tenn  goal  for  production  of  the  isomer  would  be  on 
the  order  of  100  grams  per  year  (25),  and  with  1.3  GJ  stored  per  gram  this  equates  to  130 
GJ  stored  per  year. 

In  the  case  of  most  isomers  the  energy  absorbed  during  fonnation  is  reradiated, 
but  in  rare  situations,  such  as  with  this  hafnium  isomer,  “large  differences  between  the 
spins  of  the  isomeric  and  ground  states,  (and)  differences  between  their  projections  on  the 
symmetry  axis  (K  quantum  number)”  (33 : 167)  inhibit  this  decay.  In  the  case  of  178Hf“2, 
the  spontaneous  decay  is  restricted  for  decades.  The  challenge  then  is  to  get  the  isomer  to 
release  its  stored  energy  on  demand;  a  feat  that  some  scientists  are  now  convinced  is 
achievable. 

The  Air  Force  Office  of  Scientific  Research  (AFOSR),  Defense  Advanced 
Research  Projects  Agency  (DARPA),  Department  of  Energy  (DOE),  and  the  European 
Office  of  Aerospace  Research  and  Development  (EOARD)  have  all  been  primary  sources 
of  funding  for  research  of  the  triggering  and  development  of  this  very  special  isomer  and 
a  few  others.  In  1999,  under  AFOSR  sponsorship,  an  article  entitled  “Accelerated 
Emission  of  Gamma  Rays  from  the  31-yr  Isomer  of  178Hfn2  Induced  by  X-Ray 
Irradiation,”  was  published  in  the  Physics  Review  Letters  by  C.B.  Collins  et  al.  (9)  from 
the  University  of  Texas  at  Dallas.  In  the  experiments  conducted  at  the  university,  a 
dental  x-ray  machine  set  to  endpoint  energies  of  70  and  90  keV  was  used  to  irradiate 
samples  of  178Hfn2  with  the  hopes  of  triggering  the  release  of  the  total  energy  stored 
within  the  nucleus  of  the  isomer.  A  dental  x-ray  machine  produces  bremsstrahlung 
radiation;  which  is  a  continuous  spectrum  of  x-ray  photons  generated  due  to  the 
deceleration  of  charged  particles.  In  this  case  electrons  were  the  charged  particles  being 
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decelerated  and  the  x-ray  machine  was  operated  at  15  inA.  Photon  intensities  on  the 
order  of  1010  photons/keV-s  (9)  were  produced  on  the  1  cm  diameter  target  containing  6.3 
x  1014  isomeric  nuclei.  The  consortium  (9;  10;  28 ;  29)  reported  4  ±  2 %  increases  in  the 
intensity  of  the  495  keV  emission  line  emitted  during  spontaneous  decay  due  to 
irradiation  with  photons  below  90  keV  in  energy.  The  range  of  interest  for  incident 
photon  energies  was  shortly  after  identified  to  be  between  20  and  60  keV.  This  paper 
recommended  that  further  research  be  conducted  to  obtain  better  measurements  of  the 
energy  required  to  trigger  the  decay  of  the  isomer,  and  subsequent  experiments  have 
since  tried  to  better  identify  this  requirement  as  well  as  improve  the  measurement 
techniques  that  are  used  to  identify  successful  triggering  of  the  isomer. 

The  release  of  this  initial  report,  however,  sparked  a  hot  debate  within  the 
scientific  community  partially  because  the  size  of  the  cross  section  for  photon  interaction 
reported  was  much  larger  than  allowed  according  to  theory,  and  the  4  ±  2  %  increases  in 
spectral  intensity  were  not  large  enough  to  convince  many  that  they  were  witnessing  a 
successful  triggering  event  (15).  The  reported  integrated  cross  section  (oint)  for  the 

triggering  was  10"21  cm2-keV,  which  some  scientists  claim  is  beyond  the  realm  of 
physical  possibility  (13).  The  integrated  cross  section  can  be  calculated  by  multiplying 
the  cross  section  for  photon  interaction  by  the  reaction  branch  energy  width  that  excites 
the  isomer  to  the  desired  K- mixing  level  resulting  in  decay  past  the  isomeric  to  the 
ground  state  (10).  In  turn,  this  integrated  cross  section  is  related  to  the  fractional 
enhancement  of  the  decay  rate  (f),  the  irradiating  photon  flux  ( <T> ),  and  the  lifetime  of  the 
isomeric  state  (x  =1.4  x  109  s)  as  shown  in  Eq.  [1]. 
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Under  DOE  sponsorship  a  consortium  from  Lawrence  Livennore,  Los  Alamos, 
and  Argonne  National  Laboratories  conducted  an  experiment  in  attempt  to  verify  Collins 
results  (14).  An  Advanced  Photon  Source  (APS)  was  used,  and  the  paper  released  in 
2001  (17)  reported  that  no  accelerated  emission  was  observed.  The  APS  used  produced  a 
“white”  beam  of  incident  radiation  with  intensities  ranging  from  1011  to  10 15 
photons/keV-s  across  the  20  to  60  keV  region  of  interest.  In  this  first  experiment  the 
team  at  Argonne,  with  photon  of  energy  between  20  to  60  keV  incident  on  0.03  cm2 
targets  containing  7.3  x  1014,  3.0  x  1015,  and  6.4  x  1015  isomeric  nuclei,  calculated  a  cross 
section  for  triggering  that  was  below  2  x  10"27  cm2-keV.  This  was  many  orders  below  the 
10"21  cm2-keV  results  achieved  in  the  Texas  experiments,  and  these  researchers  claimed 
that  the  results  witnessed  were  consistent  with  what  is  predicted  by  nuclear  physics  using 
an  experimental  setup  that  was  orders  of  magnitude  more  sensitive  than  the  initial  Texas 
experiments  (14). 

The  consortium  lead  by  Collins,  however,  still  had  confidence  in  the  results  of 
their  past  experiment  and  followed  it  up  with  further  experiments  using  a  monochromatic 
photon  source  (Spring- 8  in  Japan).  This  type  of  source  is  capable  of  producing  a  narrow- 
band  of  incident  x-rays  on  a  target,  as  the  energy  bandwidth  is  only  0.5  eV  wide.  The 
results  of  this  testing  indicated,  once  again,  successful  triggering  of  the  178Hfn2  and 
served  as  a  verification  of  the  results  obtained  in  earlier  testing.  In  Collins  summary  of 
the  work  accomplished  to  date  in  2001  (10),  he  reported  a  4  keV  (9  to  13  keV)  bandwidth 
for  the  absorption  of  trigger  photons  and  a  cross  section  not  much  reduced  from  the  10"21 
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cm2-keV  value  obtained  earlier,  showing  promise  of  a  self- sustained  reaction  within  a 
material  after  the  initial  triggering  event.  In  addition  to  this,  a  new  gamma-ray  line  of 
129.4  keV  was  detected  coincident  with  the  ground  state  band  213.4  keV  transition, 
signifying  a  recognizable  affect  of  photon  irradiation.  A  second  experiment  by  the  team 
at  Argonne  with  the  intent  of  studying  incident  energies  below  20  keV  did  not  detect  this 
new  line  and  produced  a  slightly  larger  cross  section  10~26  cm2-keV  than  previous 
experiments,  but  still  yielded  no  evidence  of  successful  triggering  (12). 

The  scientific  methods  utilized  by  the  researchers  in  these  studies  will  not  be 
examined  in  this  work,  but  the  work  done  by  Collins  et.  al.  (9;  10;  28;  29)  will  be  a  basis 
for  moving  on  and  addressing  issues  with  the  application  of  this  technology.  Collins  has 
gone  further  in  stating  that  the  ionization  of  L-shell  electrons  of  l78H(m2  has  a  0.2% 
chance  of  triggering  the  release  of  the  energy  stored  in  the  nucleus  (10),  and  linked 
triggered  decay  of  1 78 HP”2  with  the  Nuclear  Excitation  by  Energy  Transfer  (NEET) 
process.  This  process  begins  with  the  photo  ionization  event,  and  when  the  vacancy  is 
filled  by  the  transition  of  an  electron  from  a  higher  orbital  shell,  energy  is  transferred  to 
the  nucleus  instead  of  to  an  emitted  photon  The  benefit  of  relating  this  process  to 
triggering  is  that  the  cross  section  is  increased  due  to  the  incident  photon  being  able  to 
interact  with  the  entire  atom  instead  of  just  the  nucleus  to  initiate  decay.  Collins  reports  a 
photoionization  cross  section  at  the  Li  edge  of  7.5  x  10  20  cm2  and  links  90%  of  the 
triggering  events  to  x-ray  energies  corresponding  to  the  L-shcll  edges  for  photoionization 
of  the  electrons  of  1 7,8  Hlm2  (10).  Investigation  of  the  actual  mechanism  by  which  the 
decay  process  proceeds  will  be  avoided  in  this  study,  and  the  assumption  that  incident  x- 
rays  can  trigger  the  release  of  the  stored  energy  in  1 7xHfl2  will  be  made.  The  diagram 
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below  shows  the  natural  decay  process  beginning  with  the  emission  of  either  a  12.7  keV 
or  309.5  keV  photon.  The  actual  ranges  of  trigger  photon  energy  still  under  investigation 
are  from  9  to  12  keV,  ~40  keV,  and  ~60  keV  (34). 


Figure  1.  Energy  Level  Diagram  of  178Hfn2  Decay  (13) 

Experiments  funded  by  AFOSR  and  DARPA  are  ongoing  at  the  University  of 
Texas  at  Dallas  and  Youngstown  State,  with  aims  at  gaining  a  better  understanding  yet  of 
the  stimulated  decay,  and  work  under  EOARD  sponsorship  is  continuing  internationally 
(35).  DARPA  is  also  financing  the  production  of  this  isomer  at  several  industrial  firms. 

2.2.  Radioactive  Decay. 

2.2.1.  Forms  of  Radioactive  Decay.  There  is  a  significant  difference  between  the 
decay  products  generated  by  178Hfrl2  and  those  generated  by  conventional  fission 
reactors.  Fission  is  a  process  where  a  heavy  nucleus  is  split  by  a  fast  moving  neutron 
resulting  in  the  formation  of  new  nuclei,  neutrons,  electromagnetic  radiation,  electrons, 
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and  massless  uncharged  particles  called  antineutrinos.  The  particles  that  are  released  are 
moving  very  fast  and  their  high  kinetic  energy  values  are  the  result  of  the  release  of 
nuclear  binding  energy.  Once  these  particles  interact  and  collide  with  other  atoms  in  the 
fuel  material,  reactor  core,  or  coolant,  their  energy  can  be  transferred  to  thennal  energy. 
The  recoverable  energy  from  a  iissio  n  reaction  is  approximately  200  Me V  per  nucleus 
which  is  available  in  the  form  of  heat.  This  energy  is  distributed  among  the  fission 
fragments  and  particles  as  displayed  in  Table  3. 


able  3.  Representative  Distribution  of  Fission  Energy  (3:471) 


Energy  Source 

Fission  Energy  (MeV) 

Fission  fragments 

168 

Neutrons 

5 

Prompt  gamma- rays 

7 

Delayed  radiation 

Beta  particles 

8 

Gamma- rays 

7 

Radiative  capture  gammas 

5 

Total 

200 

As  a  comparison  only  2.446  MeV  is  given  off  from  178Hfn2  as  prompt  gamma-rays  during 
decay,  necessitating  many  times  the  number  energy  producing  events  taking  place,  when 
compared  against  fission,  if  the  energy  available  to  heat  the  propellant  is  to  be  similar. 

Fonns  of  radiation  usually  have  associated  with  them  a  Linear  Energy  Transfer 
(LET)  value  which  is  related  to  the  level  of  interaction  that  radioactive  particles  have 
with  the  medium  they  are  passing  through.  Neutrons  and  alpha  ( a  )-particles  (helium 
nucleus)  fall  into  the  category  of  high  LET  which  means  they  do  not  travel  far  before 
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interacting,  and  in  the  case  of  alpha -particles  their  energy  is  deposited  very  near  the  place 
in  which  they  are  generated.  Neutrons,  which  are  essential  for  maintaining  a  fission 
chain  reaction,  have  a  neutral  charge  which  allows  them  to  easily  penetrate  the  electron 
clouds  surrounding  other  nuclei  and  cause  other  fission  events.  On  the  other  hand,  beta 
(  [j  )-particles  (electrons  of  nuclear  origin)  and  gamma-rays  fall  into  the  category  of  low 
LET  radiation  since  they  are  not  direct  forms  of  ionizing  radiation,  and  gamma-rays  in 
particular  travel  the  farthest  through  materials  before  depositing  their  energy.  Gamma- 
rays  are  a  form  of  high  energy  electromagnetic  radiation  that  originates  in  the  nucleus  and 
is  released  in  the  fonn  of  photons  (discrete  bundles  of  energy).  They  damage  systems  by 
ionizing  atoms,  yet  they  must  first  interacting  with  anatom  which  occurs  less  frequently 
due  to  the  fact  that  they  carry  no  charge.  Charged  particles  leave  a  trail  of  excitation  and 
ionization  through  the  medium  they  transverse,  which  is  great  for  generating  heat.  In  a 
fission  reaction  the  fission  products  generated  can  be  unstable.  Decay  heat  can  be 
significant  problem  in  fission  reactors  after  the  reaction  is  stopped  since  these  unstable 
nuclides  will  decay  at  a  delayed  rate  maintaining  the  presence  of  nuclear  reactions  with 
the  core  and  perpetuating  high  core  temperatures  for  long  periods  of  time  after  shutdown 
(36:414).  Once  the  photon-particle  chain  reaction  is  halted  in  a  triggered  isomer  design, 
the  heat  production  will  stop  shortly  thereafter  simplifying  this  problem  (37)  since  the 
fonnation  of  charged  particles  will  drop  off  rapidly  with  the  cessation  of  gamma 
production 

Table  4  displays  the  quality  factor  (Q),  and  LET  value  for  the  various  forms  of 
radiation.  These  values  are  means  by  which  radiation  doses  can  be  compared.  For 
instance,  a  dosage  of  gamma- radiation  must  be  10  times  the  dosage  of  alpha -radiation 
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received  to  have  a  similar  biological  effect.  As  these  values  grow,  the  level  of  biological 
damage  grows  too  showing  that  x-rays  and  gaimna-rays  are  the  least  damaging  form  of 
radiation.  They  will  also  travel  the  farthest  from  the  source,  which  is  not  advantageous 
for  this  application 


Table  4.  Quality  Factors  for  Various  Types  of  Radiation  (36:473) 


Type  of  Radiation 

Q 

LET,  keV/micron 

x-rays  and  gaimna- 
rays 

l 

3.5  or  less 

P  -rays,  E™*  <  0.03- 
MeV 

1.7 

~7 

Naturally  occurring 
a  -particles 

10 

53 

Neutrons: 

Thennal  to  1- 
keV 

2 

7 

1-MeV 

11 

~53 

7-MeV 

7 

~23 

Radiation  dose  is  measured  in  units  of  rads,  which  corresponds  to  the  absorption 
of  100  erg/g  of  a  substance  for  any  fonn  of  radiation,  and  when  multiplied  by  Q,  the 
equivalent- dose,  measured  in  units  of  rem,  can  be  obtained  (same  is  true  of  dose  rate  and 
equivalent- dose  rate).  While  doses  on  electronic  hardware  should  be  kept  as  low  as 
possible,  it  is  proposed  that  dose  values  in  the  range  of  0.01  to  2  rad/s  (38)  should  be 
allowable. 
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Table  5.  Typical  Exposure  Numbers  (3:493) 


Exposure 

Dose-Equivalent  or  Dose- 
Equivalent  Rate 

Natural  radioactive  material  in  bones 

0.034  rem/yr 

Flight  in  Aircraft 

0.001  rem/hr  (9  km  altitude) 

Chest  x-ray  (lung  dose) 

0.01  rem 

90- day  space  station  mission 

16  rem  (NASA  estimate) 

Properly  shielding  nuclear  space 
engine 

10  rem/yr  (3.17  x  10_/  rad/s  for 
gamma-  expo  sure) 

Table  6.  Acute  Radiation  Effects  from  Whole-Body  Exposure  to  Gamma-Radiation 
_ (3:494) _ 


Acute  Irradiation 
Level  (rem) 

Acute  Somatic  Effect 

15-25 

Subtle  reduction  in  white-blood-cell  counts;  not  generally 
apparent  from  exposure  for  one  person  unless  a  blood 
sample  was  taken  before  the  exposure 

50 

Reduction  in  white-blood-cell  count  after  exposure;  the 
count  returns  to  nonnal  in  a  few  weeks 

75 

10%  chance  of  nausea 

100 

10%  chance  of  temporary  hair  loss 

200 

90%  chance  of  radiation  sickness;  moderate  depression 
of  white-blood- cell  fractions 

400-500 

50%  chance  of  death  within  30  days  without  extensive 
medical  treatment 

>600 

Lethal  to  most  people  in  3  to  30  days;  even  with 
extensive  medical  treatment,  death  is  likely  within  a  few 
months  from  infection  and  hemorrhage 

>10,000 

Lethal  within  24  hours  from  damage  to  central  nervous 
system 

2.2.2.  Attenuation,  Absorption,  and  Emission.  The  attenuation  of  x-rays  and  gamma- 
rays  in  various  materials  is  very  important  to  this  study  for  a  number  of  reasons.  In  Fig. 
2,  the  fundamental  differences  between  all  the  known  forms  of  electromagnetic  radiation 
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are  shown.  A  shorter  wavelength  corresponds  to  higher  frequency  and  higher  energy 
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Figure  2.  Electromagnetic  Spectrum  (39) 


For  a  chain  reaction  to  be  maintained  within  the  core,  trigger  photons,  in  a 
sufficient  quantity,  will  need  to  be  produced  from  interactions  between  the  released 
gamma- radiation  and  core  materials.  In  the  interest  of  shielding,  which  will  be  covered 
in  Sec.  2.2.3,  it  is  important  to  understand  the  photon-particle  interactions  taking  place  in 
order  to  have  safe,  effective  operation  Of  prime  importance,  is  the  energy  deposited  by 
this  electromagnetic  radiation  within  the  core.  This  is  the  means  of  heat  production  that 
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will  drive  propellant  temperatures  up  and  lead  to  high  levels  of  performance.  For  a 
fission  reactor  approximately  82%  of  energy  released  will  contribute  to  local  increase  in 
temperature  (3:469)  due  to  Coulombic  interactions  that  take  place  between  fission 
products  and  other  charged  nuclei.  Gamma-rays  are  not  charged  particles  and  as  a  result 
they  do  not  experience  these  Coulombic  forces  (36:90-101;  40:170-198).  Once  a 
gamma-ray  interacts  and  generates  a  charged  particle  then  its  energy  is  deposited  local  to 
the  area  of  this  interaction. 

There  are  three  manners  in  which  gamma-rays  interact  with  materials  that  must  be 
examined  in  nuclear  core  design  problems.  They  are  the  photoelectric-effect,  Compton 
scattering,  and  pair  production.  Pair  production  requires  higher  photon  energies  than 
generated  in  this  case,  so  the  focus  becomes  the  other  two  forms  of  interaction  At 
gamma  energies  between  0.3  and  10  MeV,  Compton  scattering  is  the  dominant  mode  of 
electromagnetic  interaction  Below  0.3  MeV  the  photoelectric-effect  dominates  (3:472). 
In  Compton  scattering  the  incident  photon  is  deflected  from  its  path  by  an  orbital 
electron,  conserving  both  energy  and  momentum.  The  photon  continues  on  with  less 
energy  in  some  new  direction  free  to  interact  again,  and  the  impacted  electron  recoils 
acquiring  some  kinetic  energy  from  the  photon.  During  this  process  the  only  energy 
deposited  by  the  gamma-ray  is  the  kinetic  energy  imparted  to  the  electron.  In  the  case  of 
the  photoelectric -effect,  for  which  the  probability  of  occurrence  will  increase  as  photon 
energies  decrease,  the  photon  will  be  absorbed  and  an  orbital  electron  will  be  ejected 
from  the  atom.  Thus  a  charged  particle  is  now  released,  with  energy  equal  to  that  shown 
in  Eq.  [4]. 


E  =  E  -  E 

particle  photon  ionization 


[4] 
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All  three  processes  by  which  gamma-radiation  interacts  with  materials  transfer 
energy  to  charged  particles.  Due  to  the  Coulombic  interactions  discussed  earlier,  these 
particles  moves  only  a  short  distance  before  their  energy  converts  to  heat  (36:472).  In 
order  to  capture  the  energy  of  released  gamma-rays,  the  material  and  configuration 
surrounding  the  isomer  must  be  capable  of  causing  enough  photon-particle  interactions 
that  energy  is  deposited  locally.  High  efficiencies  of  energy  deposition  will  no  doubt 
lead  to  lesser  requirements  for  power  production  by  the  fuel  materials.  The  requirement 
then  becomes  effectively  removing  this  deposited  energy  with  propellants  to  maintain 
steady  core  operating  temperatures.  In  addition,  there  is  a  need  to  maintain  a  chain 
reaction  within  the  core  which  means  trigger  photons  will  need  to  be  generated  during 
this  process  to  avoid  including  a  source  for  this  type  of  radiation  in  the  rocket  design. 

The  vacancy  left  by  the  departing  electron  from  the  photoelectric -effect  can  result 
in  the  emission  of  radiation  in  the  form  of  characteristic  x-rays  or  the  ejection  of  Auger 
electrons.  The  energy  of  these  released  x-rays  is  detennined  by  the  difference  in  binding 
energy  of  the  ejected  electron  and  the  electron  that  takes  its  place  by  transitioning  down 
from  a  higher  shell  (40 :43).  This  process  is  known  as  x-ray  fluorescence.  In  addition  to 
finding  materials  that  fluoresce  at  the  desired  photon  energy  (9  to  13,  -40,  ~60  keV),  the 
number  of  x-ray  photons  generated  is  extremely  important  to  maintaining  a  chain  reaction 
if  the  probability  of  triggering  decay  is  only  0.2%  (10),  as  stated  by  Collins.  In  addition, 
this  process  competes  with  the  ejection  of  Auger  electrons,  which  is  more  common  in 
elements  with  low  atomic  numbers  (40:45).  Auger  electrons  are  emitted  instead  of  a 
photon,  and  in  some  materials  Auger  cascades  can  occur  releasing  several  electrons  as 
transitions  continuously  occur.  Should  an  element  with  low  atomic  number  be  necessary 
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for  supplying  the  required  x-ray  fluorescence  controlling  this  process  will  require  more 
attention. 

Important  to  all  of  this  are  the  mass  attenuation  ( p,/  p  ),  mass  energy-absorption 
( \xA  /  p  ),  and  mass  energy-transfer  ( \itr  /  p  )  coefficients  of  materials  when  exposed  to 
electromagnetic  radiation  of  various  energies.  Tables  19  and  20  in  Appendix  B  contain 
some  of  these  values  for  the  materials  of  interest  in  this  study.  These  values  simply  need 
to  be  multiplied  by  the  material  density  to  yield  the  linear  coefficients  (  p  and  p_().  The 
linear  mass  attenuation  coefficient  is  a  measure  of  the  photons  statistical  interaction 
probability  per  distance  traveled,  where  1  /  p  represents  the  average  distance  traveled  by 
a  photon  before  interacting  (40).  The  linear  mass  energy- absorption  coefficient  is  a 
measure  of  the  photons  statistically  probability  of  depositing  its  energy  per  distance 
traveled.  The  linear  mass  energy- transfer  coefficient  does  not  include  the  energy  emitted 
in  the  form  of  bremsstrahlung  radiation,  and  the  relation  between  it  and  pa  is  found  in  the 
following  equation. 

[5] 

In  Turner  (40:193),  g  is  defined  as  “the  average  fraction  of  initial  kinetic  energy 
transferred  to  electrons  that  is  subsequently  emitted  as  bremsstrahlung.” 

2.2.3.  Shielding  &  Reflecting.  An  important  point  is  that  short  wavelength 
electromagnetic  wave  emission  does  not  result  in  the  release  atomic  particles  or  residual 
unstable  radioactive  states.  The  advantages  of  this  will  likely  be  evident  when 
detennining  the  amounts  of  shielding  needed  for  safe  operation.  It  is  likely  that  all  of  the 
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released  photons  will  not  be  absorbed  within  the  core,  and  some  may  make  their  way 
beyond  the  core  containment  vessel  Therefore,  a  means  of  shielding  will  be  needed  to 
protect  surrounding  equipment  from  the  adverse  effects  of  interaction  with  the  decay 
products.  As  discussed  previously,  x-rays  and  gamma-rays  easily  penetrate  most 
materials  which  will  no  doubt  be  one  of  the  major  challenges  in  the  overall  design  of  a 
Triggered  Isomer  Core  (TIC).  The  mean  free  path  for  interaction  of  a  600  keV  gamma 
photon  in  lead  is  approximately  0.7  cm,  while  the  mean  free  path  for  absorption  extends 
over  1  cm.  X-rays  penetrate  significantly  less,  as  an  approximate  mean  free  path  for  the 
absorption  of  a  10  keV  photon  in  lead  is  about  705  pm.  Mean  free  path  calculations  for 
various  materials  and  various  photon  energies  are  tabulated  in  Table  21  in  Appendix  B. 

In  a  conventional  fission  reactor  the  reaction  is  maintained  at  a  critical  state 
through  moderators,  control  rods  and  reflectors  that  maintain  the  necessary  neutron 
population  and  energies  within  the  reactor.  Preliminary  experiments  (9;  10;  28;  29) 
indicate  that  low-energy  x-rays  are  the  key  to  unlocking  the  energy  stored  in  the  178Hfn2 
isomer  so  it  is  likely  that  their  population  within  the  reactor  will  become  essential  to 
establishing  and  maintaining  a  chain  reaction.  While  x-rays  can  be  reflected  to  some 
degree  if  their  incident  angle  is  very  small  (mirrors  are  positioned  nearly  parallel  to  the 
source),  gamma-rays  are  absorbed  in  all  materials  and  no  means  of  reflecting  them  seems 
to  exist  at  present  (41).  Gamma-ray  telescopes  must  rely  on  techniques  to  simply  count 
the  gamma-rays  incident  from  a  source  in  outer  space.  The  inability  to  reflect  gamma- 
rays  becomes  important  when  the  realization  is  made  that  any  photons  not  initially 
absorbed  within  the  core  will  need  to  be  attenuated  in  the  shielding  material.  In  addition 
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there  will  be  an  energy  loss  associated  with  this  that  will  be  an  important  consideration  in 
reactor  design. 

2.2.4.  Environmental  Laws  and  Regulations .  In  designing  a  system  that  results  in  the 
release  of  radioactive  products,  engineers  will  be  faced  with  a  long,  detailed,  and  drawn 
out  process  for  the  approval  to  design,  test,  and  operate  such  a  system.  The  benefits  of 
such  an  engine  must  outweigh  the  potential  risks.  The  risks  of  flying  such  a  rocket  within 
the  atmosphere,  whether  it  is  operating  or  simply  being  transported  to  power  later  stages 
must  be  assessed  with  an  emphasis  on  environmental  impact  should  an  accident  occur.  If 
the  leap  is  made  to  manned  systems;  the  crew  must  be  adequately  protected  during 
operation  and  in  the  event  of  a  malfunction.  The  hazards  associated  with  construction, 
storage,  and  launch  pad  operation  must  be  assessed  and  a  number  of  people  present 
during  safety  reviews  and  approval  processes  must  be  convinced  that  the  design  can  be 
operated  safely.  An  example  of  the  Nuclear  Safety  Review  and  Launch  Approval 
Process  is  shown  in  Fig.  3,  and  until  this  process  is  complete  normal  launch  operation 
cannot  proceed. 


Figure  3 :  Nuclear  Safety  Review  and  Launch  Approval  Process  (42) 
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The  largest  area  of  contention  is  likely  to  be  the  classification  of  this 


technology.  While  a  nuclear  reaction  is  not  taking  place,  radiation  is  being  emitted, 
meaning  that  a  similar  process  if  not  the  same  will  be  probable  for  this  technology. 
Gamma- radiation  also  travels  farther  than  any  of  the  other  types  of  radiation  created 
during  the  fission  reaction.  Addressed  in  this  work  will  be  the  shielding  necessary  for 
utilizing  such  a  design  and  the  likely  exposure  received  by  the  rocket’s  surroundings. 

2.3.  Nuclear  Rocket  Development 

With  the  advent  of  the  first  controllable  fission  chain  reactor  in  1942,  the  far¬ 
fetched  idea  of  utilizing  atomic  energy  as  a  means  of  power  in  space  became  more  than 
just  a  dream  Throughout  history  many  programs  in  the  U.S.,  Russia,  and  European 
nations  have  made  greats  strides  in  making  this  dream  a  reality.  While  much  research  has 
been  accomplished  to  date,  nuclear  thennal  power  has  still  not  made  a  successful 
transition  to  spaceflight.  It  has  had  to  battle  strict  environmental  laws,  lack  of  funding, 
the  end  of  the  space  race,  and  support  that  disappeared  just  as  systems  reached  a  level  of 
maturity  (1;  43).  Some  of  the  most  recognizable  gas-cooled  nuclear  space  reactor 
programs  deemed  probable  for  adaptation  to  usage  of  triggered  isomer  energy  were 
researched  and  are  described  in  the  next  section 

2.3.1.  Types  of  Space  Nuclear  Reactors.  There  are  many  differences  between  the 
terrestrial  reactors  used  to  provide  electricity  here  on  earth,  the  reactors  used  to  power 
some  naval  vessels,  and  those  reactors  that  will  be  most  applicable  to  space  operations. 
The  operating  environment  and  requirements  placed  on  reactors  launched  into  space  will 
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dictate  that  they  be  compact,  capable  of  safe  modes  during  launch,  capable  of  starting  in 
space,  be  extremely  reliable,  and  make  use  of  gases  or  liquid  metals  for  coolants.  With 
the  choice  of  constructing  a  solid,  liquid,  or  gaseous-core  reactor,  the  solid-core  is  the 
simplest  and  most  well  understood  reactor  making  it  the  first  choice  for  a  space  nuclear 
power  system.  It  has  the  disadvantage  of  lower  propellant  temperature  when  compared  to 
the  other  two  designs  but,  nonetheless,  was  chosen  on  the  basis  of  documented 
engineering  experience  (3:452-459;  4:12-19;  5;  6;  18;  19;  20;  21).  With  the  type  of 
reactor  chosen  the  means  for  harnessing  its  power  becomes  the  central  area  of  focus.  The 
two  most  popular  means  by  which  this  is  done  are  known  as  nuclear  thennal  and  nuclear 
electric  propulsion.  Nuclear  thennal  propulsion  involves  heating  a  propellant  directly 
and  accelerating  the  flow  through  a  nozzle  to  produce  thrust.  Nuclear  electric  produces 
heat  energy  as  well,  but  uses  power  conversion  systems  to  create  electricity  which  can 
power  both  onboard  systems  and  propulsion  devices  such  as  ion  thrusters.  Much  greater 
levels  of  specific  impulse  are  available  with  nuclear  electric  propulsion  over  nuclear 
thennal,  which  makes  it  an  attractive  choice  for  long  duration  space  mission.  Nuclear 
thermal  propulsion,  though,  has  the  benefit  of  high  thrust  levels  and  the  ability  to  affect 
the  vehicles  velocity  on  short  time-scales,  which  make  it  a  more  attractive  choice  for 
operating  in  the  presence  of  large  gravitational  forces  and  on  quick  turn-around  missions. 
In  Sec.  2.3  the  history  of  space  nuclear  fission  reactors  is  covered  from  which  designs  for 
a  new  system  utilizing  isomer  decay  as  the  heat  source  can  no  doubt  benefit. 
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2.3.2.  Nuclear  Engine  for  Rocket  Vehicle  Application  (NERVA).  The 

ROVER/NERVA  program  is  likely  the  United  States’  most  highly  publicized  nuclear 
rocket  program.  The  ROVER  research  portion  of  the  effort  began  in  1955  as  a  program 
aimed  at  powering  Intercontinental  Ballistic  Missiles  (ISBMs)  with  the  goal  of 
developing  a  flight- rated  thennodynamic  nuclear  rocket  engine  with  75,000  pounds  of 
thrust  (44).  The  program  produced  several  reactors  that  made  it  all  the  way  through  the 
testing  stages  of  development.  The  reactors  used  uranium-235  with  graphite  moderators, 
beryllium  reflectors,  and  hydrogen  as  the  propellant  (3:452;  5;  6;  18;  44;  45).  The  fuel 
elements  were  hexagonally  shaped  with  the  uranium  fuel  dispersed  in  a  variety  of 
manners  within  a  graphite  matrix.  Coolant  channels  flowed  through  these  elements  and 
utilized  a  protective  coating  to  prevent  erosion  of  the  fuel  elements  in  order  to  maintain  a 
fixed  configuration  throughout  the  reactors  operation.  Rotary  drums  containing  neutron 
absorbing  material  were  used  to  control  the  fission  rate  while  a  surrounding  reflector  and 
pressure  vessel  maintained  the  presence  of  the  fission  products  and  necessary  reactor 
pressure.  Of  these  the  KIWI  series  of  reactors  successfully  demonstrated  the  basics  of 
building  a  nuclear  rocket.  The  program  was  able  to  overcome  the  initial  difficulties  of 
graphite  erosion  through  the  use  of  niobium  carbide  coatings,  make  a  successful 
transition  in  fuel  material  from  uranium  oxide  (UO2)  to  uranium  carbide  (UC2),  and 
eliminate  the  internal  vibrations  causing  the  fuel  elements  to  fracture.  Much  larger 
reactors  were  developed  in  the  Phoebus  series  with  the  final  test  operating  at  over  4,000 
MWt  (5).  The  coating  protecting  the  graphite  from  H2  corrosion  was  further  improved  in 
this  test  series  to  a  niobium  carbide- molybdenum  mixture.  With  the  Pewee  series  of 
reactors  aimed  and  building  a  more  compact  reactor  power  densities  were  increased  and 
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the  coating  was  improved  to  a  zirconium  carbide  material.  NASA  took  over  the  program 
in  1958  and  by  1960,  in  coordination  with  the  Atomic  Energy  Commission  (AEC);  the 
NERVA  portion  was  commissioned  along  with  the  Space  Nuclear  Propulsion  Office  to 
manage  it  (44).  This  program  went  beyond  the  current  efforts  of  researching  reactor 
design  and  was  aimed  at  producing  an  actual  space  nuclear  reactor  that  could  be  depbyed 
on  missions.  The  Nuclear  Reactor  Experiment  (NRX)  began  testing  in  1964  and  sought 
to  demonstrate  that  a  reactor  could  be  built  that  was  capable  of  withstanding  vibration 
and  shock  loadings  consistent  with  launch  vehicle  operations.  The  program  culminated 
with  the  testing  of  the  XE- Prime  reactor  which  was  the  first  to  be  fired  like  an  actual 
rocket  with  its  exhaust  end  facing  downward.  Reliability,  restart  capability, 
predictability,  controllability  and  structural  integrity  were  all  of  interest  in  this  testing  and 
the  final  reactor  design  made  significant  progress  towards  all  of  these  goals  (5).  The  XE- 
Prime  testing  achieved  power  levels  as  high  as  1,100  MWt,  but  once  again  was  faced 
with  the  problems  of  cracking  fuel  bundles  and  eroding  graphite  despite  the  evolution  of 
coating  materials  (5).  The  program  ended  in  1972  but  not  before  a  number  of  different 
reactor  designs  were  tested  at  AEC’s  Nuclear  Rocket  Development  Station  in  Jackass 
Flats,  Nevada.  The  ROVER  &  NERVA  program  generated  a  proven  reactor  design  that 
was  actually  slated  for  a  mission  to  Mars  after  the  U.S.  successfully  landed  on  the  moon. 
The  end  of  this  program  came  about  not  because  of  engineering  failures  or 
insurmountable  scientific  hurdles,  but  because  of  changes  in  govermnent  spending  and  a 
shift  in  space  priorities.  With  a  renewed  interest  in  NTR  technology,  improvements  to 
the  old  NERVA  designs  have  been  realized  with  advances  in  cooling,  material,  and 
manufacturing  technologies  (18). 
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2.3.3.  Particle  Bed  Reactor  (PBR).  In  the  early  1980s  the  idea  of  using  a  PBR  to 
increase  the  surface  area  available  for  heating  a  propellant  began  to  grow  in  popularity 
(21).  A  PBR  consists  of  a  bed  of  spherical  fuel  particles  protected  from  the  propellant 
flow  by  layers  of  graphite  and  zirconium  carbide.  The  particles  are  housed  between  two 
concentric  porous  cylinders  called  frits.  The  propellant  flows  radially  in  through  the 
outer  (cold)  frit,  is  heated  by  the  fuel  particles  in  the  packed  bed,  flows  out  radially 
through  the  inner  (hot)  frit,  then  flows  axially  out  of  the  fuel  element,  and  is  then 
expanded  through  the  nozzle  to  produce  thrust.  A  variable  number  of  fuel  elements  can 
be  included  in  a  hexagonal  arrangement,  surrounded  by  a  moderating  material  to  control 
the  reaction.  A  reflector  and  pressure  vessel  surround  this  arrangement  of  fuel  elements. 
In  1987  the  Strategic  Defense  Initiative  Office  initiated  the  Space  Nuclear  Thennal 
Propulsion  (SNTP)  program  which  assumed  control  of  the  development  of  this 
technology  (21).  This  program  was  transitioned  to  the  Air  Force  in  1991  and  only  lasted 
until  1993.  No  full-scale  engine  tests  were  conducted  as  with  the  ROVER/NERVA  tests 
series,  but  the  concept  was  proven  through  a  series  of  smaller  tests.  Fuel  tests  with 
hydrogen  reached  temperatures  as  high  as  3,000  K,  where  the  highest  temperature 
reached  in  the  ROVER/NERVA  experiments  was  2,650  K  (3 :453).  Other  achievements 
of  this  program  were  tests  of  single  fuel  elements,  criticality  experiments  for  a  prototype 
1,000  MW  core,  tests  with  power  densities  as  high  as  40  MW/liter,  various  mission 
designs,  and  verified  computer  codes  (3 :454). 

2.3.4.  CERMET.  Fast  fission  nuclear  reactors  have  also  been  given  some  attention  for 
possible  space  nuclear  reactor  designs.  These  types  of  reactors  do  not  utilize  a  moderator 
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to  control  the  fission  rate  because  they  operate  on  a  fast  fissioning  spectrum  at  energy 
levels  above  1  MeV  (3 :455)  and  rely  on  higher  uranium  enrichment  within  the  fuel 
elements  to  allow  self-  sustained  fission.  The  CERMET  reactors,  as  they  are  known, 
contain  uranium  dioxide  fuel  particles  embedded  in  a  metal  matrix  of  molybdenum  and 
uranium  or  tungsten,  rhenium,  and  uranium.  The  tenn  cermet  is  simply  used  to  classify 
such  a  metal  matrix.  Both  tungsten  and  rhenium  have  high  melting  points  and  high 
atomic  numbers  meaning  that  they  are  good  materials  for  withstanding  core  temperatures 
and  for  absorbing  the  radiation  given  off  by  the  fuel  source.  Molybdenum  has  a 
relatively  high  melting  temperature  and  can  be  used  in  the  cooler  regions  of  the  reactor  to 
decrease  the  reactor  mass.  From  1961  to  1967  Argonne  National  Laboratory  conducted  a 
program  to  develop  just  such  a  reactor  (5).  Unlike  the  ROVER/NERVA  series  no 
engines  were  built  or  tested  in  this  program  but  a  longer  operating  life,  restart  ability,  low 
sensitivity  to  temperature  cycling,  and  better  compatibility  between  hot  hydrogen  and  the 
fuel  were  all  realized  as  potential  benefits  of  the  development  of  this  technology.  Fast 
reactors  are  also  generally  smaller  and  lighter  than  thennal  reactors,  but  this  can  vary 
based  on  the  composition  of  the  cennet  matrix. 

2.3.5.  Isotope  Thermal  Thrusters .  As  a  result  of  the  work  being  done  under  the 
ROVER  program,  a  series  of  miniature  direct  cycle  nuclear  systems  were  spawned 
(46:103).  A  very  simple  form  of  heating  a  propellant  was  realized  by  utilizing  natural 
alpha,  beta,  or  gamma-radioactive  decay.  In  general  these  systems  consist  of  a 
radioisotope  capsule  with  some  type  of  decay  particle  absorbing  structure  which  is  heated 
by  the  radioactive  decay  and  able  to  transfer  its  heat  to  the  propellant  passing  over  it 


31 


before  the  propellant  is  expanded  out  a  nozzle.  The  Marquardt  Corporation  worked  on  a 
gamma-heated  system  under  USAF  sponsorship  which  utilized  thick  tungsten  heating 
elements  on  the  outside  of  the  capsule  to  capture  the  radiated  energy  (46:101).  The 
POODLE  thruster  which  had  a  and  |j  sources  also  had  both  active  and  passive  thennal 
control  systems.  The  active  control  system  which  was  able  to  expose  the  heat  source  to 
space  when  propellant  flow  was  cut  off  had  calculated  thennal  efficiencies  near  90% 
(46:102).  Though  these  were  very  low  thrust  systems  due  to  their  dependence  on  a 
natural  rate  of  radioactive  decay,  they  reinforce  the  idea  that  gamma- radiation  can  be 
captured  within  a  core  and  used  to  heat  a  propellant. 


32 


3. 


Methods  and  Theory 


The  basic  goals  when  designing  a  heat  source  for  space  propulsion,  whether  it 
powers  a  launch  vehicle  or  a  deep  space  probe,  can  be  very  similar.  Regardless  of  the 
application,  engineers  strive  to  design  a  system  that  can  effectively  and  efficiently 
transfer  heat  to  a  propellant  through  radiation  and  convection  Effective  shielding  is 
incorporated  to  allow  for  a  safe  operating  environment.  The  mass  of  this  heat  producing 
core  and  its  required  shielding  are  minimized  to  preserve  payload  capabilities.  In  the 
case  of  a  nuclear  or  triggered  isomer  core,  the  design  configuration  must  be  worked  out 
with  great  care  so  that  a  chain  reaction  can  be  maintained  and  controlled  in  such  a  fashion 
that  the  melting  temperatures  of  the  core  materials  are  not  exceeded.  Also  a  means  for 
starting  and,  in  some  instances,  stopping  the  chain  reaction  on  demand  could  be  very 
beneficial 

3.1.  Assumptions. 

3.1.1.  Isomer  Decay  Process.  Due  to  the  infancy  of  controlling  isomer  decay  and  the 
lack  of  existence  of  similar  systems,  two  general  assumptions  needed  to  be  made  in  this 
study  to  limit  the  research  in  some  manner. 

1 .  A  chain  reaction  and  the  heat  generation  rate  within  the  core  can  be  controlled  and 
steady-state  operation  is  achievable.  This  includes  a  means  for  generating  and 
maintaining  the  necessary  trigger  photon  population. 
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2.  Production  of  the  isomer  1 78Hf“2  can  be  escalated  to  at  least  the  100  g/yr  goal 

(25). 

3.1.2.  Rocket  Configuration.  There  are  a  multitude  of  systems  on  a  rocket  that  could 
be  affected  by  a  change  to  the  propulsion  system,  and  the  goal  of  this  study  was  not  to 
design  an  entirely  new  rocket  configuration  Instead,  the  current  system  was  assumed  to 
remain,  for  the  most  part,  unchanged.  Both  stages  of  the  Delta  IV-H  are  cryogenically 
fueled  and  small  modification  to  the  storage  and  handling  system,  such  as  eliminating  the 
liquid-02  segment  and  expanding  the  liquid- Hi  storage  area,  were  imagined.  In  both  the 
gas- generator  (stage  1)  and  split- expander  (stage  2)  propellant  feed  systems  used  in  the 
current  system,  regenerative  cooling  of  the  thrust  chamber  structure  pre-heats  the 
propellant  (3 :201).  For  the  purposes  of  this  study  all  the  propellant  was  assumed  to  be  in 
gaseous  state  by  the  time  it  reaches  the  turbine.  This  may  lead  to  lower  chamber 
pressure,  (3 :200)  below  what  is  desired  in  the  actual  design,  but  optimization  of  this  value 
will  not  be  addressed  in  this  study.  The  nozzle  segments,  designed  specifically  for  the 
vehicles  flight  profile,  were  not  altered  and  the  expansion  ratio  used  on  the  current 
vehicle  was  carried  over  to  this  study. 

3.2.  Nuclear  Fission  Reactor  Performance. 

Like  a  fission  reactor,  the  Triggered  Isomer  Core  (TIC)  will  benefit  from  the 
ability  to  use  a  single  propellant  for  the  dual  purpose  of  cooling  the  reacting  materials  and 
providing  the  thrust  necessary  to  propel  the  vehicle  through  space.  The  first  step  in  this 
process  was  to  use  the  infonnation  gathered  on  the  NERVA,  CERMET,  and  PBR  designs 
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(Table  1)  to  baseline  the  potential  for  improvements  in  performance  by  instituting  fission 
reactors  as  the  propulsion  sources  in  space  launch  and  orbit  maneuvering  vehicles.  The 
perfonnance  of  a  TIC  measured  against  these  fission  reactors  was  then  addressed.  Each 
reactor  was  examined,  through  use  of  the  process  outlined  in  Sec.  3.3,  to  verify  that  it 
could  satisfy  the  requirements  of  the  missions  (Table  2)  and  to  determine  what  level  of 
improvements  are  possible  over  the  current  systems.  The  data  on  the  reactor  designs 
does  not  necessarily  represent  the  best  that  can  be  achieved  with  these  designs,  but  it  was 
used  to  show  that  significant  improvements  will  still  result,  lending  weight  to  the 
argument  for  not  relying  solely  on  chemical  propulsion  for  high  thrust  missions. 


3.3.  Rocket  Fundamentals. 


3.3.1.  Nozzle  Analysis.  An  excellent  preliminary  design  process  is  presented  in  the  text 
by  Humble,  Henry,  and  Larson  (3).  The  process  outlined  there  was  used  in  this  research 
(Sections  3.3. 1-3. 3. 4)  to  perform  an  analysis  of  a  rockets  perfonnance  and  power 
requirements  knowing  the  conditions  listed  below. 


a.  Propellant 

b.  Thrust  (F) 

c.  Nozzle  expansion  ratio  ( e  ) 


f  Velocity  change  requirement  ( Av  ) 

g.  Gas  temperature  at  core  inlet  & 
outlet  ( Tx  and  T2 ) 

h.  Core  chamber  pressure  ( Pc  or  P: ) 


d.  Initial  vehicle  mass  (m(.) 

e.  Payload  mass  ( m ,  ) 


i.  Estimated  pressure  drop  through 
reactor  ( Pdrop  ) 

j .  Mach  number  at  core  outlet  ( M2 ) 
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1 :  Core 
Inlet 


2:  Core 
Inlet 


3:  Nozzle 
Exit 


Figure  4.  Position  within  Engine  as  denoted  by  Subscripts 


The  following  analysis  is  essentially  just  an  examination  of  rocket  nozzle  perfonnance 
with  the  following  assumptions. 

1 .  Isentropic  flow  (reversible  &  adiabatic)  through  the  nozzle. 

2.  Negligible  friction,  fluid  viscosity,  and  heat  transfer  to  nozzle  walls. 

3.  One  dimensional,  steady,  frozen  flow. 

4.  Propellants  are  treated  as  perfect  gases 

Propellant  specific  heats  (  c  )  were  calculated  from  equations  obtained  from  Humble, 

Henry,  and  Larson  (3 :460)  and  a  Chemkin  data  file  (47).  Knowing  the  specific  heat 
value,  the  following  expression  can  be  used  to  obtain  the  ratio  of  specific  heats  ( y  ). 


[6] 


As  stated  earlier,  the  core  inlet  temperature  in  every  case  was  assumed  to  be  300 
K  and  under  this  assumption  the  propellants  used  in  this  study  are  all  gases  at  that 
temperature.  Knowing  the  temperature  and  pressure  conditions,  and  the  flow  Mach 
number  at  the  core  outlet,  the  stagnation  properties  can  be  obtained  using  the  following 
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isentropic  relations.  The  importance  of  Mach  number  is  discussed  in  Sec.  3.3.2  of  this 
work. 


1  + 


=  P(1  + 
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Acoustic  ( aQ )  and  characteristic  (c  *)  velocities  are  now  needed  and  can  be  detennined 
via  the  following  equations. 
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The  Mach  number  at  the  nozzle  exit  must  be  iteratively  calculated  using  the  following 
equation. 


7+1 

0  =M3£  — )(1  +  ^— ^-)]2y“2 
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Now  the  ratio  of  exit  pressure  ( Pe )  to  chamber  pressure  (  Pc  )  can  be  found,  and  knowing 
either  the  chamber  or  exit  pressure  will  detennine  the  value  for  the  other. 
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The  nozzle  exit  temperature  ( Te )  is  available  through  another  isentropic  relationship  once 
the  exit  pressure  is  determined. 
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And  now  the  important  perfonnance  characteristics  such  as  specific  impulse  (  /  , ),  mass 


flow  rate  ( m  ),  and  nozzle  exhaust  velocity  ( ve )  can  be  calculated. 
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[15] 

[16] 


Pressure  thrust  has  been  neglected  in  this  ana  lysis  since  that  level  of  exactness  was  not 
needed  for  the  comparisons  made  here. 


3.3.2.  Conservation  of  Mass.  Flow  through  a  core,  as  long  as  the  flow  channels  can  be 
modeled  as  straight  circular  ducts  and  the  flow  is  steady,  can  be  analyzed  using  the 
conservation  of  mass  law.  The  importance  of  this  stems  from  the  desire  to  have  a  low 
flow  velocity  throughout  the  core  (definitely  subsonic)  and  in  the  region  of  mach  number 
0.2  to  promote  heat  transfer  and  avoid  the  loss  effects  of  flow  induced  vibration  (3 :463). 
Thus  in  this  study  the  Mach  number  of  the  propellant  exiting  the  core  was  set  at  0.2  to 
satisfy  this  requirement  leading  to  flow  at  the  inlet  in  the  range  of  0.05  to  0.16  for  the 
various  reactor  designs.  At  the  core  outlet  the  velocity  of  the  flow  ( v, )  is  needed  first. 

v2  =  M2a2  [17] 
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For  a  constant  area  duct  we  then  have  a  conserved  quantity  (  G  )  based  on  the  mass  flow 
in  equaling  the  mass  flow  out. 

G=p2v2  [18] 

The  equation  above  does  not  apply  to  the  PBR  design  since  the  flow  area  changes  along 
the  length  of  the  flow  passage,  but  with  the  documented  value  of  the  pressure  drop,  no 
major  issues  are  raised  in  the  analysis.  It  is  just  not  possible  to  calculate  Vj  and  Ml  for 
the  PBR  design  without  more  information  on  the  flow  geometry.  The  reactors  are 
assumed  to  have  a  10%  (Enhanced  NERVA),  or  53.7%  (CERMET)  pressure  drop  (  Pdrop) 
due  to  flow  losses  (3:463).  In  this  analysis  this  pressure  drop  is  given  as  a  percentage  of 
the  achieved  chamber  pressure  at  the  core  outlet  ( P2 ).  Knowing  the  pressure  at  the  core 
outlet  ( P2 )  allows  us  to  work  backwards  and  to  calculate  the  required  pressure  at  the  core 
inlet  ( P{ ). 

P^P2+PdroA  [19] 

The  same  process  shown  in  Sec.  3.3.1  was  then  used  to  find  y , ,  a, ,  P!  (at  the  core  inlet), 
and  the  flow  velocity  and  Mach  number  at  the  heater  inlet  can  be  determined  in  the  same 
manner  as  shown  in  Eqs.  [17]  and  [18].  The  same  isentropic  relations  used  in  Sec.  3.3.1 
are  also  used  to  calculate  the  stagnation  properties  of  the  propellants  at  the  reactor  inlet. 

3.3.3.  Power  Calculations.  Boundaries  have  been  set  on  the  propellant  temperature  at 
the  core  inlet  and  outlet.  The  temperature  of  the  propellant  entering  the  core  is  300  K  and 
at  the  reactor  outlet  it  is  going  to  vary  based  on  the  properties  of  the  propellant  and  the 
core’s  ability  to  transfer  heat  by  convection  to  the  propellant.  Solid  cores  are  typically 
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limited  to  the  melting  temperatures  of  the  materials  from  which  they  are  constructed.  For 
the  ceramics  used  in  the  reactors  considered  here  that  is  somewhere  in  the  range  of  3,100 
to  4,200  K  (Table  18  in  Appendix  A).  A  derivation  of  the  energy  equation  leads  to  the 
following  expression  for  detennining  the  power  necessary  to  heat  a  particular  mass  flow 
rate  of  propellant  to  the  desired  temperature.  This  is  the  power  ( Pcore)  that  must  be 

produced  and  transferred  to  the  propellant  within  the  core  for  this  technology  to  be 
successful. 

pco„=^rr(K  +  \cpdT)  P0] 

MM 

1\ 

The  above  expression  will  yield  Pcore  in  kW  and  must  be  multiplied  by  1,000  to  express 
the  value  in  watts. 


3.3.4.  System  Sizing.  The  level  of  thrust  required  is  determined  by  the  mission,  which 
will  directly  drive  the  mass  flow  rate  required  through  the  core  (Eq.  [15]).  To  increase  a 
rocket’s  thrust,  one  must  increase  the  mass  flow  rate  of  propellant  through  the  reactor.  A 
new  relation  can  now  be  used  to  show  the  link  between  the  chamber  pressure  (  Pc  )  and 


the  nozzle  throat  area  dimension  (  At  ),  which  will  directly  affect  the  engine  size  and 
weight. 


c  *  in 

p,m( 


[21] 


The  characteristic  velocity  ( c  *)  is  a  function  of  the  propellant  and  its  heated  temperature, 
and  the  use  of  this  equation  requires  the  assumption  of  isentropic,  1-D,  steady  flow  of  a 
perfect  gas.  For  a  given  expansion  ratio,  a  larger  throat  area  will  dictate  a  larger  nozzle 


40 


exit  area  and  overall  increase  in  the  size  of  the  rocket  nozzle.  As  the  rocket  nozzle  grows 
in  size  it  grows  in  weight,  so  maintaining  large  values  of  P  can  lead  to  smaller  rocket 
nozzles.  On  the  other  hand,  larger  values  of  P  also  lead  to  increases  in  the  mass  of  the 

pressure  vessel  surrounding  the  core  so  there  is  a  trade  off.  This  vessel’s  purpose  is  to 
maintain  high  pressure  throughout  the  core  and  its  mass  (m  )  can  be  estimated  with  the 

following  relation  (3 :272). 


2^(10)X 

m  = - 

9.814V 


[22] 


This  relation  is  derived  in  Humble,  Henry,  &  Larson  (3)  for  estimating  the  storage  tank 
mass  since  the  shape  of  a  reactor  pressure  vessel  is  similar  to  a  propellant  tank.  For 
nuclear  fission  reactors  a  typical  value  for  Pc  is  between  3  and  10  MPa  (3:502).  As 
shown  later,  flight  within  the  atmosphere  will  not  be  covered  in  great  detail,  but  it  is 
important  to  note  that  chamber  pressure  has  another  important  effect  on  vehicles  flown  in 
regions  of  high  atmospheric  pressure.  Flow  separation  within  the  nozzle  is  typically 
controlled  by  adjusting  the  chamber  pressure  and  nozzle  expansion  ratio  to  achieve 
values  an  exit  pressure  between  15  and  45  kPa  (3:205).  In  that  case  more  emphasis 
should  be  placed  on  this  value. 


3.3.5.  Mission  Feasibility.  The  goal  here  is  to  detennine  if  the  rocket  design  is  feasible 
for  the  intended  segment  of  the  mission.  Each  segment  has  its  own  velocity  change  ( Av  ) 
requirement,  and  each  rocket  segment  will  have  restrictions  on  propellant  mass  (mp)  and 
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gross  mass  at  liftoff  (mj).  The  amount  of  propellant  (nip)  needed  to  achieve  a  required 
Av  with  a  certain  system  is  obtainable  through  the  following  two  relationships. 


i  f  =  m,(e  ) 


mp  =  mi  -  m  f 


The  first  relationship  shown  above  is  a  variation  of  the  ideal  rocket  equation  and  it  is  used 
with  the  effects  of  the  pressure  thrust  being  neglected.  The  variable  ny  is  the  final  mass 
of  the  system  upon  achieving  the  desired  change  in  velocity  ( Av  ).  Specifically  for  this 
study  the  propellant  mass  savings  (mextra)  achieved  by  using  a  propulsion  source  capable 
of  higher  Isp  values  must  not  be  offset  by  the  mass  of  that  source  plus  the  shielding 
requirements  based  on  its  radiation  output.  The  symbol  mprop  designates  the  mass  of 
propellant  carried  by  the  current  system  and  (mextra)  is  simply  obtained  by  the  following 
relationship. 


™  extra  =  ™  prop  ~  ™  p 


The  quantity  must  be  positive  in  order  to  continue  else  the  replacement  of  a  chemical 
system  doesn’t  make  sense;  better  performance  is  desired.  The  inert  mass  fraction  ( fnert ) 
can  be  obtained  in  the  following  fashion  by  first  detennining  the  structural  mass  of  the 
rocket  (ms). 


ms  =  m  f  -  m, 


mp  +  ms 
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Another  process  generates  an  equation  for  propellant  mass  using  a  combination  of  mass 
equations  and  the  ideal  rocket  equation.  The  result  is  an  equation  with  the  following  tenn 
as  the  denominator. 

Av 

1-  finer^  >  0  [28] 

A  propellant  mass  that  is  infinite  or  less  than  zero  is  physically  impossible,  so  the 
equation  above  serves  as  a  check  for  feasibility. 

Thrust-to- weight  ( F/W )  ratio  is  also  an  important  rocket  parameter,  especially  for 
launch  systems.  A  launch  vehicle’s  initial  thrust-to-weight  ratio  must  be  greater  than  one 
for  the  vehicle  to  leave  the  ground.  Typical  values  for  this  are  in  the  range  of  1.2  to  1.5 
for  launch  vehicles,  and  usually  greater  than  0.2  for  upper  stages  (3:17-18).  If  rocket 
payloads  can  withstand  the  g- forces  of  high  thrust,  then  a  spacecraft  will  experience  a 
lesser  values  of  Av  loss  due  to  gravitational  forces.  For  the  purposes  of  this  study  a  high 
F/W  ratio  will  be  considered  a  positive  characteristic. 


3.3.6.  Rayleigh  Line  Analysis.  An  investigation  of  the  probability  for  thennal  choking 
was  also  conducted.  The  equation  used  in  this  analysis  (48 :246)  requires  the  assumption 
of  a  perfect  gas  and  constant  specific  heats,  so  an  average  value  of  the  specific  heat 
throughout  the  reactor  was  used.  First  the  stagnation  temperature  at  Mach  1  must  be 
calculated. 


T*  = 


[29] 
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The  quantity  T02/T0*  has  tabulated  values  for  given  Mach  number  (48:669).  Equation 
[30]  solves  for  the  maximum  amount  of  energy  {qmax)  that  can  be  added  to  a  fluid  before 
thennal  choking  occurs.  That  is  the  condition  that  describes  the  propellant  being 
accelerated  to  Mach  1  by  the  heating  process  alone,  and  this  will  have  serious 
repercussions  on  the  flow  through  the  nozzle.  Supersonic  flow  entering  the  nozzle’s 
diffuser  section  will  decelerate  resulting  in  large  drop  in  propellant  exhaust  velocity.  As 
the  exhaust  velocity  decreases  so  does  the  specific  impulse. 

<7imx  =TolcPavg(^-V  [30] 

1  o\ 

To  get  the  maximum  power  ( Pmax )  multiply  qmax  by  the  mass  flow  rate. 

Pmax  =  9 max  '»  [31] 


3.4.  Hafnium  Requirement. 

The  detennination  of  the  minimum  amount  of  hafnium  isomer  that  must  decay  to 
supply  the  power  required  is  calculated  with  the  following  series  of  equations.  This 
amount  of  isomer  may  not  be  sufficient  to  sustain  a  chain  reaction,  but  nevertheless 
provides  a  starting  point  for  estimating  fuel  requirements.  The  bum  rate  (cr)  is  first 
determined  knowing  both  the  power  requirements  of  the  core  ( PCOre )  and  the  energy 
stored  by  the  hafnium  isomer  (1.3  GJ/g). 


cr  = 


1.3  — *1000 


g 


[32] 


Burn  time  (th)  is  calculated  knowing  the  propellant  used  and  its  mass  flow  rate. 
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[33] 


m , 


m 


From  these  two  quantities  the  mass  (in  hj)  and  volume  (V Hf)  of  hafnium  required  can  be 
obtained  through  use  of  the  following  two  equations. 


m 


—  t  * 


Hf  ~  lb 


t.  *  cr 


[34] 


V  =  • 

y  Hf 


m 


Hf 


'  Hf 


[35] 


3.5.  Energy  Deposition  and  Fluorescence. 

Due  to  the  current  state  of  triggered  isomer  research  it  is  not  likely  that  models 
exist  for  simulating  a  chain  reaction  within  such  a  material,  and  no  models  were 
discovered  in  this  study.  Instead  the  approach  taken  in  this  study  was  to  examine  various 
materials  (those  used  in  the  fission  designs  and  other  likely  candidates),  and  assess  their 
abilities  to,  not  only,  absorb  photon  energy,  but  also  to  fluoresce  at  the  desired  photon 
energies.  Fluorescence  is  a  process  by  which  an  atom  absorbs  a  photon  and  re-emits  a 
photon  of  different  wavelength. 

For  each  of  the  nuclear  core  designs  the  power  required  by  the  core  (. Preq )  was 
assumed  to  be  the  power  required  to  heat  a  particular  mass  flow  rate  of  propellant  to  a 
desired  temperature  based  plus  power  escaping  due  to  radiation  from  the  core.  Preq  is 
actually  greater  than  what  is  needed  to  heat  the  propellant  to  the  desired  temperature  and 
the  efficiencies  for  this  process  of  these  nuclear  cores  was  estimated  from  the  data 
available  (Table  1)  by  comparing  their  power  values  against  the  power  needed  to  heat  Fk 


45 


to  the  listed  temperature  achieving  the  listed  thrust  value.  The  actual  calculation  of  these 
efficiencies  is  displayed  in  Appendix  H(H.l). 

Achieving  these  levels  of  efficiency  with  a  TIC,  however,  will  require  additional 
material  within  the  core  due  to  the  mean  free  paths  in  materials  for  gamma-ray  photons 
(Table  21  in  Appendix  B).  The  Nuclear  Engineering  Handbook,  by  Etherington,  (1 1 :7- 
66)  provides  a  process  by  which  this  efficiency  can  be  roughly  estimated.  To  carry  this 
out,  the  core  was  envisioned  to  be  a  point  source  gamma  emitter  producing  photons  of 
only  500  keV  (E0).  A  spherical  layer  of  material  is  placed  around  this  source  at  a 
distance  R ,  with  a  thickness  x  as  shown  in  Fig.  5. 


x 

Figure  5.  Core  Simplification  to  a  Spherical  Body 

In  order  to  calculate  the  energy  at  a  being  deposited  at  a  distance  R  from  the  source 
beyond  a  thickness  x  of  material  the  following  expression  can  be  used. 

H(R,x)  =  1 .062(10rV„£  A  (|ix)r°  (R,x)  [36] 

The  linear  energy-absorption  coefficient  ( (1 4 )  takes  into  account  the  most  probable  fonns 
of  photon-particle  interaction  (photoelectric -effect,  Compton  scattering,  and  pair 
production),  and  as  stated  earlier,  will  vary  with  the  energy  of  the  photon  and  type  of 
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material  being  considered.  Table  20  in  Appendix  B  lists  the  mass  energy- absorption 
coefficients  for  a  variety  of  materials  subject  to  various  incident  photon  energies.  The 
mass  energy- absorption  coefficient  is  simply  multiplied  by  the  material  density  (  p  )  to 
obtain  the  linear  energy- absorption  coefficient.  In  order  to  generate  values  for  //(the 
power  of  the  gamma  radiation  at  the  outer  surface  of  the  enclosing  layer  of  material)  the 

current  density  of  uncollided  photons  (T0 )  and  the  buildup  factors  for  energy  absorption 
(Ba)  are  required.  The  former  can  be  detennined  in  the  following  fashion. 

Se~ ^ 

T°(R,x)  =  -—  [37] 

4tiR 

The  source  strength  ( S)  was  obtained  by  dividing  an  estimated  core  power  by  500  keV 
per  photon  to  obtain  the  number  of  photons  being  produced  per  seconds.  This  should 
allow  for  a  worse  case  estimate  since  most  of  the  photons  produced  during  the  decay  of 
178  HP11 2  have  less  energy  than  this.  For  the  relaxation  lengths  (  |ix )  listed  in  Etherington 
(11)  and  with  knowledge  of  the  linear  attenuation  factors  for  various  materials  (Table  18), 
layer  thickness  for  the  listed  relaxation  lengths  can  be  found.  Each  thickness  will  have 
associated  with  it  a  buildup  factor  ( Ba ).  This  process  is  demonstrated  in  Appendix  I  for 
lead  and  tungsten.  Since  data  on  tungsten  was  not  available  in  the  tables  it  was  obtained 
by  iterating  between  two  of  the  material  listed  in  the  table.  The  fact  that  buildup  factors 
are  linear  functions  of  the  atomic  number  allows  for  this.  The  volume  of  material 
absorbing  the  radiation  was  then  used  to  provide  the  mass  of  material  required  and  also 
detennine  the  total  power  being  radiated  beyond  the  surrounding  layer  of  material.  For 
various  amounts  of  material  it  is  possible  to  detennine  a  probable  efficiency  for  this 
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energy  release.  Core  masses  will  likely  need  to  be  increased,  by  at  least  these  amounts 
calculated,  if  energy  deposition  within  the  core  is  to  be  achieved. 

A  materials  ability  to  absorb  incident  gamma-radiation  and  subsequently  emit  x- 
ray  photons  of  the  energy  needed  to  sustain  a  chain  reaction  within  the  TIC  was 
examined  using  only  the  listings  of  x-ray  transition  energies  for  various  materials  (49).  A 
material’s  x-ray  photon  emission  energy  (or  transition  energy)  is  determined  by  the 
difference  between  the  ionization  energy  of  ejected  electron  and  that  of  the  electron  that 
transitions  down  from  a  higher  orbit.  The  energies  listed  in  the  tables  provided  by  the 
National  Institute  of  Standards  and  Technology  (NIST)  are  not  necessarily  representative 
of  achievable  photon  energies  as  some  may  occur  with  very  low  probability.  For  this  to 
work  the  probabilities  for  ionizing  just  the  right  orbital  electron  followed  by  the  transition 
producing  photons  in  the  isomer’s  range  of  sensitivity  to  triggering  will  need  to  be  very 
high  since  only  a  limited  number  of  gamma-photons  are  emitted  from  178Hfn2  as  it 
decays.  Materials  were  chosen  from  these  listings  as  potential  sources  of  trigger  photons 
in  a  TIC  based  on  a  combination  of  their  potential  for  emitting  photon  within  the  desired 
range  and  their  level  of  attenuation.  It  would  be  desirable  to  have  the  photons  produced 
reach  the  isomeric  material  in  the  core  before  being  absorbed  near  the  point  of  their 
emission  (34). 

3.6.  Radiation  Shielding. 

The  shielding  analysis  carried  out  here  focuses  on  shielding  for  biological  and 
equipment  protection  purposes  and  not  on  heat  shielding  which  is  used  to  protect  the 
pressure  vessel  surrounding  the  core  from  excessive  heat  (36:548).  Radiation  shielding 
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in  general  is  dependent  on  many  factors  such  as  the  power  and  type  of  radiation  emitted, 
the  distance  between  the  radiation  source  and  those  components  sensitive  to  radiation, 
and  the  level  of  sensitivity  of  those  components  (3 :494).  For  the  purposes  here  the 
mission  is  unmanned  eliminating  the  need  to  protect  humans,  but  in  the  case  of  the  1 st 
stage  vehicle  shielding  may  be  necessary  to  prevent  a  harmful  release  into  the 
enviromnent.  The  assumption  of  radiation  hardened  electronics  and  payload  is  made  in 
this  study,  to  a  level  able  to  withstand  doses  up  to  0.01  rad/s  (38).  Safe  doses  for  human 
are  in  the  area  of  1  prad/s  (2:66),  and  this  level  of  radiation  exposure  will  be  examined  to 
view  how  shield  weights  can  vary  with  the  level  of  radiation  exposure  tolerable. 

3.6.1.  Nuclear  Fission  Shielding.  Shielding  from  the  harmful  effects  of  a  nuclear 
fission  reaction  is  somewhat  different  from  shielding  against  a  gamma  release  alone. 
Present  in  the  fission  reaction  are  fast  moving  neutrons  and  alpha  particles  which  are  high 
LET  forms  of  radiation  These  particles  have  high  probabilities  of  interaction  with  other 
particles  within  the  medium  they  are  passing  through.  It  is  also  not  always  desirable  to 
absorb  all  of  this  radiation  so  a  variety  of  materials  can  be  used;  one  is  beryllium  to 
reflect  neutrons  and  help  maintain  their  population  within  the  core,  and  another  is  lithium 
hydride  used  to  slow  down  neutrons  as  well  as  take  some  part  in  their  attenuation.  A 
basic  flat  shield  design,  shown  in  Fig.  6,  is  capable  of  reducing  the  gamma  ray  flux  by  a 
factor  of  0.00105  and  the  neutron  flux  by  a  factor  of4.0(10)'9  (3:495).  The  image  shown 
is  a  cross-sectional  view  showing  the  thickness  of  the  shields  material  layers  of  beryllium 
(Be),  tungsten  (W),  and  lithium  hydride  (LiFL). 
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Radiation  _ 

Be  W  LiH2 

(18  cm)  (5  cm)  (5  cm) 

Figure  6.  Typical  Fission  Reactor  Shadow  Shield  (3 :496) 

The  dimensions  of  this  shield  can  be  adjusted  to  achieve  further  levels  of  reduction,  but 
as  shown  such  a  shield  would  have  a  shield  loading  factor  of  3,500  kg/mr,  and  calculation 
of  the  shield  weight  would  simply  require  multiplying  the  shields  surface  area  by  this 
value.  For  case  where  a  shadow  shield  is  place  on  top  of  the  reactor  to  shield  the 
remainder  of  the  spacecraft  from  escaping  radiation,  the  radius  of  this  shield  can  be 
approximated  to  have  the  same  radius  of  the  reactor  it  is  shielding. 


3.6.2.  Gamma-Ray  Shielding 

The  method  used  for  calculating  the  radiation  shielding  necessary  for  protection 
from  a  gamma- ray  release  was  taken  from  Turner  (40:368)  and  Lamarsh (36:549).  This 


method  also  treats  the  core  as  a  point  source  emitter.  The  dose  rate  ( D  )  received  at  a 
distance  r  from  the  source  can  then  be  determined  knowing  the  mass  energy- absorption 
coefficient  for  air  ( ji  ,  /  p  )  and  the  rate  of  energy  being  released  in  the  form  of  gamma- 
rays  ( CE ). 


D  = 


CE  \lA 
4 nr2  p 


[38] 
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The  values  for  the  mass  energy- absorption  coefficient  vary  based  on  photon  energy,  but 
for  air  the  value  is  pretty  consistent  across  the  range  of  100  to  600  keV  photons.  An 
important  consideration  in  the  actual  design  will  be  the  medium  separating  the  spacecraft 
from  the  energy  source  since  this  will  affect  the  dose  rate  by  altering  (  p, A  /  p  ).  Once  the 
dose  at  a  distance  without  shielding  is  known  the  shield  relaxation  length  (  p,x)  can  be 
determined  with  the  following  relation  including  a  buildup  factor  ( B )  that  takes  into 
account  the  scattering  of  photons  within  the  material. 

(ur=-kjg(-A)  [39] 

DoB 

The  process  involves  detennining  both  B  and  \ix  iteratively  since  both  are  unknowns. 


Relaxation  Lengths 

Figure  7.  Buildup  Factor  vs.  Relaxation  Lengths  (40) 
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A  guessed  value  of  B  is  input  into  Eq.  [39]  and  varied  until  the  calculated  value  of  the 
relaxation  length  the  value  for  B  match  as  shown  in  Fig.  7.  Lead  was  chosen  as  the  shield 
material  in  this  study.  It  is  likely  that  the  cost  of  producing  lead  shields  would  be  less 
than  some  of  the  other  metals  considered  in  this  study,  but  it  may  not  be  the  best  option 
for  an  actual  design  based  on  its  low  melting  temperature. 

Invariably  the  most  important  value  to  detennine  is  the  shield  loading  (measured 
in  kg/m2)  which  allows  one  to  calculate  the  shield  mass  knowing  only  the  surface  area  of 
the  shield.  Shield  loading  is  obtained  by  multiplying  the  shield  thickness  by  the  density 
of  the  shielding  material.  An  important  point  to  make  here  is  that  this  analysis  is  valid 
only  for  a  monodirectional  beam  normally  incident  on  a  slab  shield  (36:553)  since  the 
buildup  factors  used  were  generated  specifically  for  such  an  application.  Shields  of 
various  different  shapes  can  use  the  same  method  but  require  different  buildup  factors. 
Consider  also  that  178Hfn2  is  not  a  source  for  a  monodirectional  beam  of  gamma-radiation 
when  viewing  these  results.  The  hope  is  that  the  worst  case  scenario  is  captured  by 
assuming  the  emission  of  600  keV  photons,  where  in  reality;  the  photons  released  will 
range  between  12.7  and  574-keV  (see  Fig.  1). 

Table  23  in  Appendix  E  shows  the  trends  established  in  this  study  of  shielding  a 
gamma-ray  source.  Quantities  such  as  the  distance  from  the  source  (r),  the  allowable 
dose  rate  ( D ),  the  distance  between  the  shield  and  the  source,  and  the  percentage  of 
power  radiated  to  the  surroundings  were  varied  for  a  constant  area  shield  to  detennine 
their  affects  on  shield  thickness,  loading,  and  ultimately  mass.  The  vales  shown  in 
Appendix  E  are  generic  and  not  specifically  associated  with  any  of  the  core  designs  in 
this  study.  These  examples  were  merely  canied  out  to  gain  an  understanding  of  shield 
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placement  and  the  variation  of  shield  weights  that  could  be  expected  in  the  actual 
designs. 

3.7.  Individual  Reactor  Design  for  Fission  and  Isomer  Cores. 

The  process  of  designing  a  TIC  will  look  first  at  the  materials  that  would  be  the 
best  candidates  to  make  up  the  core  based  on  their  history  of  use  in  such  an  application, 
material  properties,  and  material  interaction  with  electromagnetic  radiation.  Choice  of  a 
particular  core  configuration  depends  on  mission  requirements,  heat  transfer  capabilities, 
necessary  alteration,  and  manufacturability.  The  TIC  core  concept  was  applied  only  to 
the  PBR  design. 

3.7.1.  NERVA.  The  NERVA  design,  as  discussed  early,  incorporates  approximately 
300  hexagonally  shaped  fuel  elements  with  the  uranium  fuel  dispersed  in  a  variety  of 
manners  within  a  graphite  matrix  (3 :452).  The  rate  of  fission  within  the  core  is  controlled 
by  rotary  drums  in  the  radial  reflector  region  of  the  reactor  and  coolant  channels  run 
through  the  fuel  elements  to  cool  the  reactor.  An  analysis  of  replacing  the  uranium  fuel 
particles  with  178Hfm2  was  not  possible  in  this  study,  nor  recommended.  Little  detail  was 
found  regarding  the  arrangement  of  the  fuel  within  the  graphite  matrix,  and  realistic 
estimates  of  the  achievable  temperatures  in  this  type  of  design  are  not  possible  for  a  TIC 
at  this  time.  The  maximum  temperature  achieved  in  the  NERVA  test  series  was  2,361  K 
(3 :457)  since  the  configuration  used  is  not  the  one  best  suited  for  heat  transfer.  Materials 
used  in  the  NERVA  core  had  melting  temperatures  well  above  3,000  K  (Table  1 1). 
Specific  sizing  equations  were  not  available  without  an  in  depth  analysis  (3 :476)  of  the 
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fission  process,  and  that  was  deemed  unnecessary  here.  Instead,  in  the  case  of  the  fission 


design  the  core  volume  ( Vcore )  and  mass  (mcore)  were  estimated  by  using  quantities  such 
as  the  power  required  (Preq)  and  power  density  (Pden)  of  the  core. 
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3.7.2  CERMET.  The  CERMET  design  is  similar  to  the  NERVA  reactor  in  that 
hexagonal  shaped  fuel  elements  with  straight  coolant  channels  were  utilized.  Other  than 
differing  on  the  principal  of  fast  fission,  the  CERMET  design  also  utilized  uranium  fuel 
particles  embedded  in  a  metal  matrix.  In  addition  to  the  uranium  fuel  the  matrix 
materials  consisted  of  molybdenum,  tungsten,  and  rhenium.  Both  tungsten  and  rhenium 
have  melting  temperatures  on  par  with  the  carbide  compounds  used  in  many  reactor 
designs,  and  are  also  both  very  good  attenuators  of  gamma-rays.  As  with  the  NERVA 
design,  surface  area  for  heat  transfer  is  not  maximized  in  this  design  resulting  in  a  reactor 
that  would  not  live  up  to  the  true  potential  of  the  heat  generated  within  it.  For  the  fission 
reactor  design  there  are  linear,  least-squares  curve  fits  for  their  CERMET  sizing  (3:489). 


=  0.0034 Preq  +  20.79 

[42] 

=  0.0067/^+41.418 

[43] 

These  were  used  to  estimate  the  core  dimensions  for  a  fission  core.  The  core  volume 
(VCOre)  was  then  obtained  with  the  equation  Vcore  =  Rcore2*Hcore  and  the  core  density 
( P core)  can  be  used  to  calculate  the  core  mass  ( mcore )  by  the  equation  mcore  =  pcore  *VCOre  ■ 
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3.7.3.  Particle-Bed  Reactor.  The  unique  nature  of  the  PBR  design  allowed  for  a  closer 
look  at  the  heat  transfer  taking  place  within  the  core,  and  the  method  employed  (3 :491) 
starts  with  the  analysis  of  a  single  pellet.  Assuming  that  some  of  the  pellets  are  hafnium 
isomer  and  some  are  another  metal  suitable  to  absorbing  the  energy  given  off  by  the 
isomer,  and  knowing  the  maximum  propellant  temperature  achieved  in  the  core,  the 
temperature  at  the  pellet- coating  interface  can  be  detennined  making  the  following 
assumptions. 

1 .  Every  pellet  is  surrounded  by  a  protective  layer  of  zirconium  carbide  with  a 
thickness  of  1  pm. 

2.  Steady  state  heat  generation  and  removal  can  be  achieved.  This  results  in  a  steady 
state  process  where  the  energy  of  the  photons  released  during  the  decay  is 
absorbed  in  the  core  materials,  conducted  to  the  materials  outer  surface  (in  this 
case  the  pellet’s  outer  surface),  and  finally  transferred  to  the  propellant  by 
convective  heating. 

3.  The  spectrum  of  photons  released  contains  or  generates  trigger  photons  within  the 
core  will  trigger  the  decay  of  additional  atoms  of  the  isomer. 

4.  The  flow  speed  through  the  reactor  is  low  enough  to  allow  boundary  layer  effects 
at  the  pellets  outer  surface  to  be  neglected.  Subsequently  there  is  not  a  large 
temperature  difference  between  the  temperature  of  the  propellant  at  the  location 
of  maximum  temperature  within  the  core  and  the  temperature  of  the  pellet’s  outer 
surface. 
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The  equation  governing  conduction  used  in  this  process  is  derived  from  Fourier’s  law  for 
conduction  and  applied  to  the  spherical  pellets  of  the  design. 

AT 

Qx  -  ~kA—  [44] 

Ax 

For  metals  needing  protection  from  the  flow  of  fuel  the  coating  of  the  pellets  used 
in  this  study  was  assumed  to  be  made  up  entirely  of  zirconium  carbide  (ZrC),  but  in 
actuality  there  are  likely  a  couple  of  thin  layers  of  graphite  separating  the  zirconium 
carbide  from  the  inner  material  Though  important  for  actual  design  the  impact  of  this  is 
assumed  to  be  negligible  in  the  heat  transfer  analysis  conducted  here.  The  pellet  size  was 
set  by  assuming  a  pellet  inner  diameter  of  300  microns  (hafnium  isomer  or  other 
material)  and  an  outer  diameter  of  500  microns  which  includes  the  layers  of  protective 
coating.  Knowing  this  both  the  pellets  outer  surface  area  (represented  by  A  in  Eq.  [44]) 
and  the  pellet  volume  can  be  determined.  An  approximation  was  made  that  88%  (3 :491) 
of  the  volume  occupied  by  a  bed  of  pellets  is  actually  the  pellets  with  the  remaining  12% 
being  the  voids  between  pellets. 

The  process  began  with  setting  the  desired  propellant  temperature  to  2,400  K. 

This  is  below  the  melting  point  of  hafnium,  a  necessity  of  this  design  if  metallic  hafnium 
is  used,  and  far  enough  below  this  melting  point  to  allow  the  hafnium  isomer  to  exceed 
this  temperature  without  melting.  In  the  instance  that  hafnium  carbide  can  be  used  the 
propellant  temperature  is  escalated  to  3,500  K  in  this  study,  since  higher  values  would 
exceed  the  limits  of  the  specific  heat  equations  being  utilized.  From  this  the  power 
necessary  to  achieve  these  temperatures  can  be  obtained,  as  well  as  the  rocket 
perfonnance  characteristics  such  as  specific  impulse,  the  initial  mass  savings,  and  overall 
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mission  feasibility  (Sec.  3.3.5).  This  calculated  core  power  requirement  for  heating  was 
applied  to  two  hypothetical  reactor  efficiencies  (82%  and  67%),  to  simulate  the  power 
lost  from  the  core.  The  efficiency  of  82%  is  based  on  the  data  provided  for  the  PBR 
(Table  1)  and  its  calculation  is  shown  in  Appendix  H  (H.l).  The  value  of  67%  was  used 
to  show  the  affect  of  decreased  efficiency  on  core  mass.  Knowing  the  power  requirement 
(P req)  allows  for  determining  the  reactor  dimensions  using  Eqs.  [45]  and  [46]  (3:487) 
shown  below.  These  two  equations  calculate  the  radius  and  height,  respectively,  of  a  19- 
element  PBR.  Equations  for  two  other  possible  configurations  are  displayed  in  Appendix 
D.  Different  reactor  dimensions  and  masses  result  from  each  configuration,  and  a  choice 
is  made  usually  to  select  the  lightest  reactor,  taking  into  account  that  the  smallest  radial 
dimension  will  also  lead  to  a  less  shadow  shield  requirement. 

Rcore  =  -2.655(10 )"12/*  +  8.946(10)-X  -1. 1703(10)-^ 

[45] 

+  7.427  (1 0)“3  Pr2eq -2.2955Preq  +313.34 

Hcore  =  -4.027 (1 0)  ~5  P*q  +  0.  \A21Preq  + 17.9883  [46] 

Knowing  core  radius  and  height  one  simply  needs  to  apply  the  relations  shown  earlier 
(Sec.  3.7.2)  to  compute  core  volume  and  mass.  For  TIC  application  this  mass  value  does 
not  accurately  represent  the  changes  that  will  need  to  occur  in  this  design. 

The  fission  design  of  the  PBR  as  described  earlier  (see  Sec.  2.3.3)  when  applied 
as  a  TIC  design  would  no  longer  require  the  moderator  blocks  to  be  of  a  material  suitable 
for  slowing  down  neutrons  produced  during  fission.  The  material  that  is  used  for  the 
pellets  in  the  bed  should  also  be  changed.  While  the  hafnium  isomer  will  be  present  in 
the  form  of  some  of  these  pellets,  in  a  quantity  that  will  be  set  by  the  needs  of  the  system, 
the  remainder  of  these  pellets  can  be  constructed  of  another  material  possessing  the 
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desired  properties.  A  TIC  should  consider  using  materials  that  are  good  attenuators  of 
photons  and  possibly  producers  of  trigger  photons  as  a  result  of  the  gamma-radiation 
initially  absorbed.  The  original  designs  have  core  densities  around  1,600  kg/m3  (3 :490). 
Deduced  from  a  calculated  core  volume,  the  densities  of  uranium-carbide,  zirconium- 
carbide,  and  beryllium,  and  values  for  the  pellet  size  was  the  following  fractional  break 
up  of  the  core;  8%  of  the  core  by  volume  is  fuel  pellets  with  coating,  58%  of  the  core  is 
moderator  material,  and  the  remainder  is  flow  passages  for  the  propellant  and  spacing. 

Of  the  8%  that  is  pellets,  21.6%  of  this  volume  is  actual  fuel,  with  the  remaining  78.4% 
being  coating  material.  Appendix  H  (H.3)  displays  these  calculations  which  are  based  on 
the  recommendation  of  l/3rd  of  the  core  by  mass  being  fuel  pellets  and  their  coating  with 
the  remaining  two -thirds  being  moderator  (34).  With  these  values  a  mass  estimate  for  a 
core  altered  for  isomer  decay  can  be  obtained  replacing  materials  as  desired. 

The  number  of  pellets  (Np)  within  the  core  will  detennine  how  much  power  is 
available  per  pellet  in  this  analysis.  This  is  detennined  by  dividing  the  mass  of  the  fuel  in 
the  core  by  the  mass  of  an  individual  fuel  pellet. 
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The  power  absorbed  and  available  for  heating  the  propellant  was  then  calculated  by 

P 


Qx  = 


N, 


[48] 


Rearranging  Eq.  [44]  then  yields 


T  =  T  + 


QAX 

kA 


[49] 
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For  the  case  of  metallic  hafnium,  the  propellant  temperature  T  is  set  at  2,400  K  and  the 
temperature  of  the  fuel  pellets  surface  area  is  calculated.  The  value  of  the  thennal 
conductivity  was  varied  to  view  the  impact  on  energy  transfer  since  thennal  conductivity 
values  could  change  in  the  instance  that  other  materials  are  used  (Appendix  H  (H.4.2)). 
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4.  Results 


4.1  Propellants.  The  performance  requirements  for  a  rocket  engine  in  this  study 
were  set  by  the  need  to  perform  the  missions  of  the  1st  and  2nd  stages  of  the  Delta  IV-H 
(Table  2).  The  first  step  taken  in  this  process  was  to  verily  what  propellants  would  be 
suitable  to  carry  out  these  missions  in  a  monopropellant  reactor  design.  Each  nuclear 
core  (NERVA,  CERMET,  and  PBR)  was  examined  with  various  propellants  (Table  7) 
and  judged  on  the  basis  of:  its  ability  to  perform  the  mission  (Eq.  [28]),  the  mass  of 
propellant  required  for  the  mission  (Eq.  [25]),  and  the  potential  for  thennal  choking  (Eq. 
[31]). 

Table  7.  Specific  Impulse  (Isp)  Values  for  Various  Propellants 


Engine 

NERVA 

CERMET 

PBR 

Stage 

1st 

2nd 

1st 

2nd 

1st 

2nd 

Temperature 

(K) 

2,361 

2,507 

3,200 

Expansion 

Ratio 

21.5 

285 

21.5 

285 

21.5 

285 

h2 

798 

873 

824 

903 

941 

1037 

CH4 

322 

379 

332 

395 

374 

447 

CO 

216 

234 

223 

244 

252 

277 

co2 

186 

213 

192 

222 

217 

252 

n2 

216 

234 

222 

244 

251 

277 

02 

203 

222 

210 

232 

239 

265 

c3H8 

269 

265 

281 

282 

330 

331 
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In  both  cases  hydrogen  is  the  only  propellant  capable  of  perfonning  the  mission. 
This  is  no  surprise  considering  that  H2  has  the  lowest  molecular  weight.  This  leads  to 
higher  values  of  specific  impulse  and  exhaust  velocity  than  heavier  propellants.  Methane 
(CH4)  was  the  one  other  propellant  possessing  a  specific  impulse  value  near  that  of  the 
chemical  rocket;  calculated  to  be  447  s  in  a  2nd  stage  PBR  design.  There  would  be  little 
value  in  its  utilization,  since  reactor  masses  and  the  necessity  of  shielding  would  likely 
increase  the  overall  rocket  mass  without  a  reduction  in  fuel  usage.  Even  for  propellant 
temperatures  above  3,200  K,  specific  impulse  values  for  the  other  propellants  do  not 
increase  significantly,  and  there  is  no  promise,  for  an  application  such  as  this,  in  the  other 
propellants  studied  here.  Figure  8  shows  the  specific  impulse  variation  as  a  function  of 
temperature  for  a  2nd  stage  PBR  with  a  nozzle  expansion  ratio  of  285: 1.  These  values 
were  calculated  with  the  method  described  in  Sec.  3.3.1. 


Figure  8.  Isp  variation  with  temperature 
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The  performance  that  results  from  using  fb  will  be  examined  in  Sections  4.2  and 
4.3,  but  the  choice  of  H2  as  a  rocket  systems  propellant  brings  along  some  important 
engineering  considerations.  In  order  to  store  hydrogen  onboard  a  spacecraft  it  must  be 
cooled  to  its  liquid  state  at  20.39  K,  which  requires  special  attention  to  the  coolant 
storage  and  transportation  system  within  the  rocket.  On  a  thennal  rocket  this  means  that 
heat  from  the  reactive  core  must  be  prevented  from  raising  the  temperature  of  the 
propellant  prior  to  its  release  from  the  storage  tank.  In  both  stages  of  the  Delta  IV-H 
considered  here,  cryogenic  hydrogen  is  already  utilized  proving  that  storage  issues  can  be 
handled,  yet  the  possibility  that  escaping  radiation  could  prematurely  heat  the  propellant 
must  be  given  some  thought  during  shielding  design.  Also,  metallic  hafnium,  proposed 
to  be  the  source  of  energy  in  this  system,  rapidly  absorbs  hydrogen  at  temperatures 
around  973  K  to  form  the  compound  HfH  (50).  If  the  isomer  can  produced  in  the  fonn  of 
hafnium  carbide  (HfC),  this  will  cease  to  be  a  concern. 

There  may  be  potential  in  some  other  light  propellants  not  considered  here,  but  for 
the  purpose  of  addressing  the  design  of  a  TIC  and  making  comparisons  against  fission 
reactors,  the  choice  of  hydrogen  is  made.  Further  time  could  be  invested  in  the 
performance  and  properties  of  hydrogen  at  the  operating  temperatures  achieved,  but  this 
should  not  be  a  large  concern  during  preliminary  design  since  the  fission  reactors  tested 
in  the  past  made  use  of  this  propellant  (3 :457;  5). 

4.2  Analysis  of  Stage  1. 

The  nuclear  fission  reactors  examined  (NERVA,  PBR,  and  CERMET)  all  vary 
with  respect  to  their  abilities  to  achieve  high  propellant  temperatures  and  maintain  high 
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values  of  chamber  pressure  throughout  the  reactor.  Included  here  are  the  thennal 
efficiencies,  some  probable  values  of  the  core  densities  (3 :490)  and  power  densities.  The 
achievable  power  density  of  the  CERMET  reactor  was  not  found,  but  is  said  to  be  the 
highest  of  all  the  reactors  due  to  the  absence  of  moderator  material  (3 :489).  See 
Appendix  H  (H.l)  for  the  calculation  of  these  efficiencies. 


Table  8.  Nuclear  Core  Approximations 


Reactor 

Thermal 

Efficiency 

Core  Density 
(kg/m3) 

Power  Density 
(GW/m3) 

NERVA 

83%* 

2,300 

1.57 

CERMET 

90%* 

8,500 

>  40  (variable) 

PBR 

82%* 

1,600 

40  (variable) 

Shown  in  Table  9  are  the  performance  values  of  the  three  nuclear  fission  designs 
in  a  scenario  replacing  the  3  RS-68  engines  on  the  Delta  IV-H  with  a  nuclear  heat  source. 
As  stated  earlier  the  specifics  of  the  system  design  are  not  the  focus  of  this  study  and  the 
chamber  pressure  and  its  effect  on  system  mass  and  perfonna  nee  will  not  be  examined 
here,  beyond  noting  that  the  lower  chamber  pressures  will  lead  to  lighter  pressure  vessels 
and  larger  nozzle  dimensions  for  these  nuclear  designs  (Sec.  3.3.4).  The  values  listed  are 
indicative  of  what  can  be  achieved  with  these  reactors,  but  in  no  way  are  meant  to 
represent  the  limitations  of  these  designs.  The  specific  impulse  values  for  each  fission 
reactor  listed  in  Table  9  are  those  generated  in  this  study,  with  the  nozzle  expansion  ratio 
of  the  Delta  IV-H  systems,  along  side  those  that  are  obtained  from  Table  1  (3 :457)  which 
are  denoted  by  the  asterisk. 
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Table  9.  1 st  Stage  Nuclear  Performance  Comparison. 


Engine 

RS-68 (22) 

NERVA 

CERMET 

PBR 

Propellant 

LOX/LH2 

lh2 

lh2 

lh2 

Chamber 
Pressure  (MPa) 

9.21 

6.89 

4.14 

6.18 

Specific  Impulse 
(s) 

420 

798  825* 

824  930* 

941  971* 

Propellant 
Required  (kg) 

598,800 

420,360 

412,380 

380,210 

m  (kg/s) 

Varies 

1,270 

1,229 

1,077 

Power  (MW) 

41,794 

43,623 

51,832 

Engine  Mass 
(kg) 

-6,604  each 
{3}  (51) 

61,227 

Engine  Volume 
(m3) 

? 

26.62 

" 

Full  Shield  Mass 
(kg) 

>144,289 

F/W 

1.4 

<1.39 

-7,314 

Propellant  requirements  for  all  of  the  nuclear  designs  drop  off  appreciably  due  to 
the  increased  values  of  specific  impulse  available.  As  a  check  of  the  specific  impulse 
values  calculated  the  specific  impulse  was  also  detennined  in  this  work  using  the 
documented  expansion  ratios  only  varied  from  the  tabulated  data  on  average  by  4%. 
Appendix  H  (H.2)  displays  the  manner  in  which  specific  impulse  can  be  obtained 
utilizing  the  procedure  described  in  Sec.  3.3.1. 

A  significant  issue  manifests  when  one  tries  to  estimate  the  size  of  a  reactor 
needed  to  provide  a  thrust  sufficient  to  launch  a  vehicle  of  this  magnitude  into  space. 
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The  sizing  estimates  based  on  reactor  power  (Sec.  3.7)  for  the  CERMET  and  PBR 
reactors  don’t  apply  for  power  requirements  above  a  few  thousand  megawatts,  and  in 
order  for  the  power  requirements  of  an  individual  reactor  to  be  this  low;  at  the  fewest  1 5 
CERMET  reactors  or  26  PBR  (of  the  37  element  configuration)  would  be  needed.  These 
numbers  were  obtained  by  dividing  the  core  powers  listed  in  Table  9  by  3,000  MW  and 
2000  MW  respectively.  This  would  likely  never  be  considered  as  an  option  regardless  of 
the  source  of  heat  energy  due  to  shear  complexity.  Therefore,  values  for  a  CERMET  and 
PBR  engine  and  shield  were  not  calculated  for  this  stage.  The  mass  and  volume  of  a 
NERVA  engine  can  be  approximated  for  this  application  via  the  method  described  in  Sec. 
3.7  and  is  displayed  in  Appendix  H  (H.5),  and  it  can  be  shown  that  to  perfonn  this 
mission  with  one  engine  it  would  have  to  be  approximately  26.6  nr  in  volume  and  weight 
61,226  kg.  This  is  assuming  that  it  operates  with  100%  efficiency  at  the  established 
power  density  of  1,570  MW/m3.  If  the  linear  trend,  provided  up  through  values  of  2,000 
MW  (3 :489),  holds  for  NERVA  reactors  of  much  greater  power  levels  then  these  values 
should  be  accurate  to  within  a  few  percent. 

Shielding  necessary  to  provide  protection  against  the  radiation  being  released  was 
based  on  the  shield  loading  factor  of  3,500  kg/m2  and  an  approximated  shield  surface 
area  of  41  m2  (approximating  the  reactor  height  to  be  twice  the  value  of  the  radius).  In 
actuality  the  shield  surface  area  would  likely  be  larger  as  this  value  is  a  probable  surface 
area  of  the  reactor  itself,  and  the  shield  loading  would  most  likely  increase  since  the 
power  levels  considered  here  are  many  times  larger  than  those  of  the  reactor  designs. 

This  proposed  shield  would  surround  the  cylindrical  reactor  concentrically  and  on  top, 
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resulting  in  a  shield  mass  greater  than  144,289  kg  offsetting  the  benefits  of  the  high 
specific  impulse  value  obtained. 

During  such  a  1st  stage  ascent,  a  rocket  using  1 78H(m2  as  the  source  of  energy 
would  require  somewhere  between  1 1  and  14  kg  of  the  isomer  to  undergo  decay,  not  to 
mention  the  additional  hafnium  needed  in  the  configuration  to  maintain  the  chain  reaction 
throughout  the  ascent.  The  equations  of  Sec.  3.4  were  used  to  obtain  these  values.  At  the 
estimated  near  tenn  rate  of  production  (100  g/yr)  it  would  take  140  years  to  produce  this 
much  hafnium  isomer.  For  this  reason  alone  it  would  not  be  prudent  to  consider  further 
the  construction  of  a  TIC  for  such  high  thrust  missions,  not  to  mention  of  the  material 
cost.  No  cost  figures  were  determined  in  this  study,  but  the  price  is  certain  to  be  high  due 
to  its  limited  availability  and  the  nature  of  the  process  through  which  it  is  obtained.  The 
method  described  in  Sec.  3.6.2  cannot  be  used  to  estimate  shield  thickness  or  shield 
loading  factor  for  power  levels  as  high  as  predicted  here.  The  level  of  radiation  released 
extends  beyond  the  data  that  supports  this  study  (Fig.  7).  Dense  materials  are  needed  for 
the  attenuation  of  gamma-rays,  and  shield  loading  factors  will  be  greater  than  those 
obtained  in  the  following  analysis  of  stage  2  which  are  in  the  region  of  1,300  to  1,400 
kg/m2.  Once  again  there  will  be  a  large  impact  on  the  weight  savings  earned  with  high 
specific  impulse  values,  and  Fig.  9  shows  that  the  choice  to  not  include  shielding  would 
not  be  a  smart  one.  A  source  operating  at  giga-watt  power  levels  and  only  radiating  5% 
of  the  energy  generated  would  result  in  exposure  rate  as  high  as  17  rad/s  at  1,609  m  (one 
mile)  from  the  source.  At  this  rate  subtle  effects  of  radiation  could  be  detected  in  people 
after  just  a  few  seconds  of  exposure  (Table  6),  and  this  would  not  be  acceptable  by 
today’s  standards  of  safety.  Dose  is  shown  in  a  log  scale  in  Fig.  9,  for  two  levels  of 
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radiation  release,  showing  how  dose  with  distance  from  the  source  fa  11s  off  according  to 


the  inverse  square  law  (Mr2). 


Figure  9.  Radiation  Dose  vs.  Distance 

4.3  Analysis  of  Stage  2. 

Estimates  of  engine  mass  and  volume  are  now  obtainable  for  all  of  the  fission 
core  designs  since  the  core  power  levels  required  for  this  mission  fall  within  the  effective 
range  of  the  equations  used  (Eqs.  [42],  [43],  [45],  and  [46]).  These  estimates  are 
displayed  in  Table  10  below.  The  calculations  are  carried  out  in  Appendix  H  (H.2,  H.5, 
and  H.6)  for  the  various  core  designs.  Fission  shielding  was  discussed  in  Sec.  3.6.1  and 
individual  core  design  in  Sec.  3.7. 


67 


For  the  2nd  stage  application  all  three  of  the  fission  reactor  designs  are  predicted  to 
result  in  a  significant  mass  savings  over  the  chemical  system  used  today.  One  major 
advantage  that  would  result  from  applying  nuclear  technology  to  a  mission  such  as  this 
would  be  the  ability  to  place  shielding  only  between  the  spacecraft  and  the  reactor  itself, 
since  radiation  sent  out  into  the  space  environment  is  of  little  concern  Of  course  this 
may  raise  issues  if  the  spacecraft  is  to  be  used  near  other  spacecraft,  satellites,  or  space 
stations,  since  shield  weights  will  escalate  as  the  surface  area  requirement  grows. 


Table  ’ 

10.2nd  Stage  N 

uclear  Perfonnance  Comparisons 

Engine 

RL10B-2 

(22) 

NERVA 

CERMET 

PBR 

(19- element) 

Propellant 

LOX/LH2 

lh2 

lh2 

lh2 

Chamber 
Pressure  (MPa) 

3.21 

6.89 

4.14 

6.18 

Specific 
Impulse  (s) 

462.4 

873 

825* 

903  930* 

1,037  971* 

Propellant 
Required  (kg) 

26,529 

17,310 

16,870 

15,162 

m  (kg/s) 

24.24 

12.84 

12.41 

10.82 

Power  (MW) 

502 

483 

629 

Engine  Mass 
(kg) 

~  277  (52) 

736 

600 

547 

Engine  Volume 
(m3) 

? 

0.32 

0.0706 

0.342 

Shadow  Shield 
Mass  (kg) 

1,511 

553 

1,303 

Vehicle  F/W 
(including 
shield) 

0.27 

0.35 

0.35 

0.37 

Engine  F/W 
(including 
shield) 

40.5 

5.0 

9.7 

6.1 

Mass  Savings 

_ (kg) _ 

7,249 

8,783 

9,794 
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The  rocket  mass  and  velocity  requirements  of  this  type  of  mission  appear  to  be  the 
area  where  a  nuclear  source  or  TIC  could  have  the  largest  impact  and  would  be  best 
suited  to  be  introduced.  The  2nd  stage  mission  has  a  much  lower  thrust  requirement 
(1 10,000  N),  which  is  less  than  the  documented  capabilities  of  the  original  reactor 
designs  (Table  1).  It  is  very  important  to  show  that  these  designs  can  achieve  the 
indicated  propellant  temperature,  at  this  leve  1  of  thrust.  According  to  the  data  provided  in 
Table  1,  the  PBR  design  is  capable  of  exhaust  velocities  of  9,526  m/s  (Eq.  [16])  and 
propellant  temperatures  of  3,200  K.  This  results  in  a  mass  flow  rate  of  35  kg/s  in  order  to 
achieve  the  documented  333,617  N  of  thrust  (Eq.  [15]).  Mass  flow  rates  for  the  PBR  on 
this  missionare  around  10.82  kg/s  (comparing  values  in  Table  1  and  Table  10),  thus  the 
detennination  is  made  the  propellant  temperatures  assumed  at  the  exit  are  indeed 
achievable. 

Environmental  concerns  are  minimized  because  the  core  could  be  kept  inactive  at 
low  altitudes,  and  the  short  nature  of  the  mission  would  guarantee  that  such  a  system 
could  be  tested  for  its  length  of  operation  in  space  on  the  ground  prior  to  space  flight.  A 
mishap  during  launch  would  also  not  generate  the  radioactive  particles  that  result  in  long 
lasting  contaminants  at  the  crash  site,  making  this  a  safer  choice  than  a  fission  reactor. 

The  minimum  quantities  of  178Hfil2  needed  to  supply  the  energy  for  this  mission 
(as  determined  by  the  method  in  Sec.  3.4)  are  also  much  less,  ranging  from  432  to  555 
grams,  for  the  various  reactor  designs.  Production  rates,  while  still  not  sufficient,  would 
not  have  to  improve  much  over  current  goals  for  this  to  be  a  near  tenn  reality.  The 
importance  of  this  is  that  if  production  goals  of  lOOg/yr  are  reached,  then  the  amounts  of 
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hafnium  isomer  needed  to  perfonn  these  types  of  missions  could  be  attainable  within  4-5 


years. 

4.4  Designing  the  TIC. 

With  the  basics  laid  for  the  level  of  improvement  that  can  be  expected  for  the  2nd 
stage  mission  by  making  the  change  to  a  fission  based  power  source,  we  can  now  begin 
to  put  together  the  specifications  for  a  TIC  and  estimate  its  level  of  perfonnance.  The 
operating  principals  for  a  TIC  would  be  much  the  same  as  for  the  fission  designs 
considered  here,  hence  their  importance  to  this  work.  The  propellant  needs  to  be  heated 
to  temperatures  as  high  as  possib  le,  such  that  the  maximum  specific  impulse  values  can 
be  obtained.  This  will  take  place  by  having  the  propellant  flow  through  a  core  of 
materials  that  serve  the  purposes  of  releasing  energy,  absorbing  energy,  providing 
protection  to  other  materials,  providing  shielding  for  the  rocket  or  surroundings,  and 
supporting  the  physical  arrangement  of  the  core  itself. 

4.4.1  Shielding  a  Gamma-Ray  Producing  Source.  Table  23  in  Appendix  E  contains 
the  data  generated  in  a  series  of  calculations,  for  generic  cases,  varying  the  important 
factors  of  shield  design.  These  factors  are  the  radiated  power,  the  distance  of  the  shield 
from  the  source,  the  distance  of  the  equipment  from  the  source,  and  the  allowable  level  of 
exposure  for  whatever  material  is  being  protected.  As  a  result  the  analysis  described  in 
Sec.  3.6.2  leads  to  the  following  conclusions  about  shield  design  for  this  type  of  source. 

A  sample  calculation  for  a  19-Element  TIC  design  at  82%  efficiency  is  displayed  in 
Appendix  H  (H.4.1). 
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1 .  The  placement  of  the  shield  with  respect  to  the  source  does  not  affect  the  level 
of  exposure  to  the  radiation.  The  shield  can  be  placed  as  geometric  constraints  of  the 
rocket  design  allow.  For  shielding  that  extends  beyond  the  simple  geometry  of  a  shadow 
shield,  such  as  a  cylindrical  design  or  spherical  design,  the  shield  should  be  as  close  to 
the  source  as  allowable  to  reduce  its  size  and  weight.  For  this  study  a  shield  distance  of 
0.5  m  was  selected. 

2.  For  a  decreasing  allowable  exposure  rate,  to  levels  as  low  as  1  prad/s,  the 
values  of  the  buildup  factor  and  relaxation  length  were  off  the  charts  used  in  this  study 
and  only  placing  great  distances  between  the  source  and  the  equipment  (out  to  50  m) 
reduced  the  values  enough  to  pennit  computation.  For  manned  and  vehic  les  bearing 
sensitive  equipment  this  could  be  a  large  issue  since  heavier  shields  will  be  required 
(reducing  the  mass  savings).  This  is  an  area  for  further  investigation. 

3.  Equipment  placed  farther  from  the  source  required  less  shielding  due  to  the 
inverse  squared  law  of  radiation  falloff  (Fig.  9). 

4.  The  increase  in  shield  weight  for  a  50%  efficient  source  as  compared  to  a  95% 
efficient  source  was  only  94  kg.  This  is  the  difference  in  the  mass  values  in  column  13 
and  column  1  of  Table  23.  The  achievement  of  lower  levels  of  core  efficiency  with 
increased  losses  to  radiation  is  not  a  large  concern  for  the  shield  design.  This  could  be 
especially  beneficial  in  the  case  that  waste  heat  needs  to  be  radiated  from  the  core,  with 
steady  requirements  for  rocket  protection. 

4.4.2  Material  Selection.  Due  to  the  conditions  within  the  core  and  the  nature  of  the 
electromagnetic  radiation  being  released  in  the  decay  process  of  178Hfn2,  careful 
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consideration  needs  to  be  given  to  the  choice  of  materials  used.  Materials  within  the 
core,  most  importantly,  must  be  resistant  to  high  temperatures  if  any  increases  in 
perfonnance  are  to  be  seen  They  must  also  maintain  good  physical,  thermal,  and 
mechanical  properties  in  a  harsh  environment  of  high  pressures  and  temperatures  for  the 
period  of  time  that  the  mission  dictates.  That  internal  environment  will  be  somewhat 
different  from  a  fission  reactor  since  the  presence  of  the  highly  ionizing  radiation  is 
absent,  namely  the  alpha  particles  and  neutrons.  Instead  of  materials  that  moderate, 
reflect,  and  absorb  neutrons,  the  TIC  will  need  materials  that  are  good  absorbers  of  short 
wavelength  electromagnetic  radiation  (gamma-rays).  For  the  purpose  of  maintaining  a 
steady-state  reaction  within  the  core,  means  for  producing  the  necessary  trigger  photons 
is  needed,  and  certain  materials  may  possess  the  ability  to  satisfy  this  need  through 
fluorescence.  Every  material  within  the  core  may  not  satisfy  all  these  requirements,  but 
the  combination  of  these  traits  between  the  materials  included  is  essential  for  the 
realization  of  this  idea.  The  background  discussion  of  the  nuclear  fission  reactors  in  Sec. 
2.3  mentions  several  materials  that  could  be  essential  to  a  TIC  design,  and  Table  1 1  lists 
those  important  properties  for  selecting  those  that  would  be  the  best  candidates  for  the 
design.  Metallic  hafnium  has  a  melting  temperature  of  2,506  K,  and  that  temperature  is  a 
limiting  factor  in  this  core  design  should  hafnium  carbide  not  be  attainable.  Regardless, 
no  materials  with  lower  melting  temperatures  than  2,506  K  were  selected  for  inclusion  in 
the  core.  The  assignment  of  a  low,  medium,  or  high  value  to  a  material’s  level  of  photon 
absorption  is  simply  a  distinguisher  for  the  materials  researched  in  this  study.  Those  with 
a  high  mark  were  the  best  attenuators  examined  in  this  study  and  so  forth.  Information 
for  some  elements/materials  could  not  be  found. 
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Table  11.  Recommendation  on  Materials 


Element 

Use  in  Fission 
Reactor 

Be 

Reflector 

C 

Support/Ceramic 

Formation 

A1 

Pressure  Vessel 

B 

Control  Drums 

Zr 

“ 

ZrC 

Protective 

Coating 

Nb 

“ 

NbC 

Protective 

Coating 

Mo 

Attenuation 

Hf 

Control  Drums 

HfC 

Ta 

“ 

W 

Attenuation 

Re 

Attenuation 

Os 

“ 

Pt 

“ 

Level  of 
photon 
absorption 

Potential 
source  of 
trigger 
photons 

Melting 

Temperature 

(K) 

Low 

? 

1,560 

Low 

? 

3,800 

Low 

No 

933 

Low 

? 

2,573 

Low 

No 

2,128 

? 

? 

3,813 

Medium- 

Low 

No 

2,750 

? 

? 

4,033 

Medium 

No 

2,896 

Medium 

Yes 

2,506 

? 

? 

4,173 

High 

Yes 

3,290 

High 

Yes 

3,683 

High 

Yes 

3,443 

High 

Yes 

3,323 

High 

Yes 

2,041 

Pb 


Medium 


Yes 
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The  best  absorbers  of  radiation  at  short  wavelengths  (x-rays  and  gamma-rays)  appear  to 
be  the  dense  transition  metals  such  as  hafnium,  tantalum,  tungsten,  rhenium,  and  the 
highly  toxic  metal  osmium.  Those  materials  in  the  shaded  blocks  are  the  most  favored 
candidates  for  a  TIC  since  they  posses  a  combination  of  the  desired  thennal  properties 
and  can  serve  either  to  absorb  radiation  or  fluoresce  as  a  result  of  it.  The  category 
entitled  “potential  source  of  trigger  photons”  needs  clarification  at  this  time.  The 
assigmnent  of  a  “yes”  or  “no”  in  this  column  simply  indicates  whether  or  not  x-ray 
transition  energies  were  found  for  this  element  in  the  range  of  interest  (9  to  13,  -40,  and 
~60  keV).  The  ability  to  maintain  a  chain  reaction  in  the  core  without  the  aide  of  an 
external  photon  source  would  be  of  great  value  to  this  system.  Materials  have  been 
identified  that  posses  the  desired  melting  temperature  and  attenuation  properties,  yet  also 
have  electron  transition  that  produce  photons  of  the  energies  viewed  to  trigger  the  decay 
of  178Hf”2.  The  probabilities  for  these  transitions  have  not  been  established  in  this  study, 
and  it  is  recognized  that  none  of  these  materials  may  actual  satisfy  the  requirements  for 
this  type  of  design. 

4.4.3  Configuration  Selection.  The  choice  was  made  in  this  study  to  apply  the 
properties  of  a  TIC  to  the  PBR  design  configuration.  The  requirements  of  the  mission 
being  that  of  an  upper  stage  vehicle  and  the  ability  of  the  particle-bed  design  to  transfer 
the  greatest  amount  of  heat  to  the  propellant  were  heavy  lactors  in  this  decision.  Large 
surface  areas  are  required  for  heating  propellants  such  as  (H2,  NH3,  and  TbO)  to  absorb 
heat  directly  from  a  reactor  (53).  It  is  also  possible,  based  on  the  individual  nature  of  the 
fuel  particles  in  this  design,  to  conduct  a  more  thorough  analysis  of  the  heat  transfer 
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taking  place  as  discussed  in  Sec.  3.7.  This  provides  a  more  realistic  view  of  the 


relationship  between  propellant  and  fuel  particle  temperatures  within  the  core.  A  detailed 
analysis  of  the  heat  transfer  taking  place  in  a  design  such  as  this  taking  into  account  local 
heat  variation,  and  radiation  is  not  needed  at  this  time.  When  the  triggering  of  nuclear 
spin  isomers  is  finally  proven  and  shown  to  be  capable  of  maintaining  a  chain  reaction 
such  an  analysis  would  definitely  be  required.  If  hafnium  carbide  is  used  this  becomes  a 
mute  point  for  this  analysis  since  the  chosen  temperature  of  3,500  K  is  far  below  the 
melting  temperature  of  that  material. 

As  mentioned  in  Sec.  3.7.3  the  thennal  conductivity  was  varied  to  identify  its 
effect  on  interior  pellet  temperature.  The  plot  below  shows  that  there  is  a  minimum 
required  thermal  conductivity  of  the  coating  material  if  the  desired  propellant 
temperatures  are  to  be  reached  without  exceeding  the  melting  temperature  of  the  fuel 
pellet.  The  data  supporting  this  plot  has  been  tabulated  in  Appendix  H  (H.4.2). 
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Figure  10.  Effect  of  Thennal  Conductivity  on  Propellant  Temperature 
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Since  the  melting  temperature  for  metallic  hafnium  is  2,506  K,  coating  material  thermal 


conductivities  should  exceed  a  value  of  1.0  W/m-K.  The  thennal  conductivity  of  1.9 
W/m-K  provided  by  Humble  (3:492)  for  zirconium  carbide  coated  pellets  was  used  in 
this  study. 

With  the  choice  of  2,400  K  (metallic  hafnium)  and  3,500  K  (hafnium  carbide)  as 
the  desired  static  propellant  temperature  at  the  core  outlet,  the  following  values  (Table 
12)  of  specific  impulse,  minimum  reactor  power,  and  mass  savings  are  achievable  for  this 
mission.  The  mass  savings  shown  does  not  include  the  subtraction  of  the  chemical 
engines  mass,  the  addition  of  the  TIC  mass,  or  the  required  shielding.  The  method  of 
Sec.  3.3.1,  used  for  the  fission  core  calculations  (Appendix  H  (H.2.1  and  H.2.2)),  was 
also  used  to  generate  the  values  listed  in  Table  12.  The  flow  geometry  of  the  PBR  is 
assumed  to  remain  unchanged  leading  to  chamber  pressures  no  different  from  the  fission 
core  designs. 


Table  12.  First  Iteration  Core  Properties 


Propellant 

Temperature 

(K) 

Isp 

(S) 

Minimum 

Core 

Power 

(MW) 

Mass 

Savings 

(kg) 

2,400 

881 

422 

9,340 

3,500 

1,090 

548 

11,960 

For  the  case  of  the  isomer  in  metallic  fonn,  realizing  that  422  MW  of  power  is  needed  to 
heat  liquid  H2  to  the  desired  temperature  providing  the  required  thrust  of  1 10,000  N;  two 
core  efficiencies  were  examined  to  view  the  effect  on  the  overall  size  of  the  core  and  the 
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shielding  requirements.  While  the  specifics  of  the  core’s  arrangement  was  not  examined 
in  this  study,  it  is  understood  that  some  of  these  pellets  will  be  constructed  of  the  hafnium 
isomer,  and  the  total  mass  of  isomer  residing  within  the  core  will  have  to  be  adjusted 
depending  on  the  amount  needed  to  sustain  a  chain  reaction  at  the  desired  rate.  Sample 


calculations  for  metallic  hafnium  pellets  are  shown  in  Appendix  H  (H.4). 


Table  13.  Mass  of  Materia’ 


is  Replacing  Pe’ 


lets  in  a  19-Element  PBR 


Core 

Efficiency 

82% 

66.7% 

Number  of 
Pellets 

392,610,596 

418,157,176 

Mass  of 
Coating  (kg) 

136 

144 

Mass  of 
Tungsten  (kg) 

104 

111 

Mass  of 
Hafnium  (kg) 

74 

79 

Mass  of 
Tantalum 
(kg) 

92 

98 

Mass  of 
Rhenium  (kg) 

117 

124 

As  stated  earlier  only  somewhere  in  the  region  of  432  to  555  grams  of  isomer  would  be 
need  to  undergo  decay  to  supply  the  energy  for  this  mission,  so  conceivably  the 
percentage  of  pellets  within  the  core  that  are  constructed  of  the  isomer  could  be  very 
small.  It  is  also  conceivable  that  the  remainder  of  the  pellets  making  up  the  core  could  be 
a  mixture  of  materials,  but  the  results  shown  in  Table  13  are  rough  estimates  of  the  total 
pellet  mass  if  a  majority  of  the  pellets  are  any  one  of  the  materials  listed.  Including 
hafnium  in  this  list  assumes  that  pure  metallic  form  of  the  element  in  a  non- excited  state 
is  being  utilized  in  addition  to  the  isomer.  The  coating  composition  was  not  varied  in  this 
study  and  its  mass  remains  the  same  regardless  of  the  material  used. 
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Table  14.  Mass  of  Materials  Replacing  Moderator  in  a  19-Element  PBR 


Core 

Efficiency 

82% 

66.7% 

Core  Volume 
(m3) 

0.3212 

0.3421 

Mass  of 
Tungsten  (kg) 

3,506 

3,734 

Mass  of 
Hafnium  (kg) 

2,480 

2,641 

Mass  of 
Tantalum  (kg) 

3,102 

3,304 

Mass  of 
Rhenium  (kg) 

3,916 

4,171 

The  replacement  of  the  moderator  within  the  core  with  materials  that  were 
identified  as  the  best  absorbers  of  high  energy  electromagnetic  radiation  and  as  possible 
sources  of  the  desired  fluorescence  will  lead  to  much  heavier  core  designs  than  tho  se 
predicted  by  using  a  typical  PBR  core  density.  The  densities  of  neutron  absorbing 
material  such  as  beryllium  and  lithium  hydride  are  much  less  than  those  of  the  metals 
shown  in  Table  14.  A  sample  calculation  is  shown  in  the  beginning  of  Appendix  H 
(H.4). 

In  the  search  for  lighter  elements  possessing  the  desired  properties  stated  above, 
no  elements  were  identified  as  potential  sources  of  the  desired  photon  emission  with 
melting  temperatures  even  approaching  2,000  K.  Provided  a  source  for  photon  emission 
can  be  found  amongst  the  heavier  elements,  the  designer  will  be  faced  with  the  challenge 
of  exposing  the  isomer  to  these  photons  before  they  are  absorbed  within  the  material  that 
produces  them  (34). 

A  shadow  shield  constructed  of  lead  and  placed  0.5  m  from  the  core,  with 
allowable  levels  of  exposure  on  the  back  side  set  at  0.01  rad/s,  will  result  in  the  shield 
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masses  shown  in  Table  15  for  two  different  levels  of  core  efficiency.  As  shown  by  the 
small  masses  calculated  below  the  ability  to  simply  use  shadow  shields  is  of  great  benefit 
a  core  designed  for  2nd  stage  applications.  Appendix  H  (H.4.1)  contains  the  calculations 
for  the  82%  efficiency  case. 


Table  15.  Shadow  Shielding  a  19-Element  TIC-PBR 


Radiated 

Power 

Core 

Radius 

(m) 

Tolerable 
Dose  Rate 
(rad/s) 

Shield 

Thickness 

(kg) 

Shield 

Loading 

(kg/m2) 

Shield 

Mass 

(kg) 

18% 

0.3558 

0.01 

0.1208 

1,377 

548 

33% 

0.3428 

0.01 

0.1263 

1,440 

535 

In  this  instance  the  shield  mass  is  actually  less  for  the  lower  efficiency  case,  but  this  is 
due  solely  to  the  fact  that  the  increase  in  core  power  to  account  for  the  radiation  losses 
has  resulted  in  a  core  with  smaller  radial  dimensions  lessening  the  shadow  shield  are. 
According  to  the  conduction  law  applied  (Eq.  [44])  and  the  assumption  of  an  equal 
distribution  of  the  energy  released  from  the  isomer  within  the  packed  bed,  the  melting 
temperature  of  metallic  hafnium,  and  of  course  tungsten,  will  not  be  exceeded  in  this 
design.  This  is  a  very  simplified  analysis;  however,  assuming  that  the  core  as  it  is 
designed  would  result  in  such  heat  flow  behavior.  In  actuality,  the  system  would  need  to 
be  modeled  at  a  very  detailed  level  in  coordination  with  the  actual  core  design. 

Presuming  that  a  chain  reaction  can  be  initiated,  much  higher  power  densities  could  be 
achieved  within  the  core,  with  effective  heat  removal  governing  the  design  (34).  For  this 
mission,  however,  a  core  power  density  similar  to  that  shown  in  Table  16  would  be 
sufficient,  and  the  amount  of  hafnium  isomer  required  to  supply  this  power  will  be  a 
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function  of  its  susceptibility  to  triggering  and  the  trigger  photons  production  within  the 


core. 


Table  16.  Actual  Core  Attributes  for  the  2,400  1 

C  Case 

Core 

Efficiency 

Total 

Core 

Power 

(MW) 

Power 

Density 

(GW/m3) 

Pellet  Core 
Temperature 
(K) 

Power 
Absorbed  / 
Pellet  (W) 

82% 

443 

1.602 

2,472 

1.075 

66.7% 

633 

1.850 

2,468 

1.009 

As  a  final  comparison,  probable  TIC  designs,  utilizing  a  configuration  similar  to 
that  of  the  82%  efficient,  19-element  PBR  design,  are  examined  next  to  provide  a 
comparison  with  the  data  tabulated  on  the  NERVA,  CERMET,  and  19-element  PBR 
fission  reactor  designs  and  RL10B-2  (Table  10).  The  TIC  as  shown  is  of  three  different 
constructions.  One  consists  of  replacing  the  pellets  and  moderator  material  of  the  fission 
design  with  all  metallic  hafnium,  and  hafnium  isomer  pellets  are  dispersed  throughout  the 
bed  as  needed.  This  would  be  the  lightest  TIC  configuration.  The  other  two  replace  the 
pellets  and  moderator  material  with  the  metal  tungsten  and  hafnium  carbide.  The  fact 
that  some  of  the  pellets  within  the  tungsten  design  will  actually  be  hafnium  should  not 
impact  the  mass  estimates  by  much  since  both  hafnium  and  hafnium  carbide  are  less 
dense  than  tungsten.  Also  note  that  the  hafnium  carbide  case  is  capable  of  higher  core 
temperatures.  All  of  the  new  core  designs  result  in  a  level  of  improvement  over  the 
baseline  system  making  a  strong  argument  for  their  place  as  space  propulsion  options. 
Provided  shadow  shields  can  support  the  mission  and  the  equipment  on  board  is  hardened 
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to  radiation  exposure  levels  of  0.01  rad/s,  very  significant  mass  savings  can  be  achieved 


with  any  design  as  shown  in  Table  17. 


Table  17.  2nd  Stage  TIC  1 

Mass  Comparisons 

Core 

Composition 

Metallic 
Hf  TIC 

Tungsten 

TIC 

Hafnium 

Carbide 

TIC 

Specific 
Impulse  (s) 

881 

881 

1,090 

Propellant 
Required  (kg) 

17,189 

17,189 

14,569 

Power  (MW) 

515 

515 

668 

Engine  Mass 
(kg) 

2,689 

3,746 

2,691 

Engine  Volume 
(m3) 

0.3212 

0.3212 

0.3342 

Shadow  Shield 
Mass  (kg) 

548 

548 

460 

Vehicle  (F/W) 

0.31 

0.31 

0.35 

Engine  F/W 

3.46 

2.61 

3.56 

Mass  Savings 
(kg) 

6,380 

5,323 

9,086 

A  few  notes  on  these  numbers.  The  core  masses  have  been  calculated  based  on  core  sizes 
predicted  for  a  fission  reaction  The  amount  needed  will  depend  heavily  on  the  level  of 
gamma  absorption  required  for  heating  and  the  mass  numbers  calculated  in  the  analysis 
shown  in  Appendix  I  for  high  levels  of  efficiency,  suggest  that  these  masses  could  be 
less.  Secondly,  the  thrust- to- weight  values  are  not  maximum  values  since  increases  in 
propellant  will  allow  for  increases  in  mass  flow  rate,  and  increases  in  thrust  for  the  same 
core  weight.  The  vehicle  thrust- to- weight  (F/W)  has  been  calculated  by  setting  the  new 
vehicle  weight  to  a  value  equal  to  the  original  vehicle  mass  minus  the  mass  savings.  The 
manner  in  which  this  savings  is  applied  is  up  to  the  user  and  the  type  of  mission  required. 
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The  choice  can  be  made  of  carrying  near  the  current  system’s  fuel  load  to  increase  the 
velocity  change  capabilities  of  the  rocket.  More  Av  can  mean  faster  transit  times,  further 
range  in  space,  longer  duration  missions,  or  increases  in  the  number  or  type  of  maneuvers 
carried  out  by  the  spacecraft.  Durability  of  spacecraft  could  also  be  improved  with  the 
additional  weight  being  used  to  design  more  rugged  components,  better  shielding  from 
the  space  environment,  and  systems  capable  of  being  used  for  longer  periods  of  time. 

Due  to  the  launch  costs  of  today,  the  most  important  benefit  of  such  a  system  would  be 
the  increase  in  payload  mass  that  can  be  put  into  space.  With  a  PBR  design  relying  on 
isomer  decay  from  a  hafnium  carbide  core  an  additional  9,086  kg  can  be  boosted  into 
space.  This  would  nearly  double  the  payload  limit  for  the  current  mission.  Satellites  tend 
to  weigh  from  a  few  hundred  to  a  few  thousand  kilograms  (22:107-1 10),  so  if  the  payload 
bay  can  accommodate  it  and  the  fuel  economy  will  allow  it  the  number  of  satellites 
placed  into  orbit  by  this  single  vehicle  could  help  reduce  the  number  of  launches  need  for 
placement  of  our  space  assets.  With  planned  missions  to  the  moon  and  requirements  for 
construction  in  space  (8),  such  as  with  the  international  space  station,  boosting  greater 
payloads  weights  into  orbit  is  of  vital  importance  in  the  near  future. 
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5.  Conclusions  and  Recommendations 


In  this  study  three  nuclear  fission  reactors  and  a  conceptual  core  dependent  on 
triggered  isomer  decay  were  examined  in  the  role  of  providing  thrust  for  the  1st  and  2nd 
stages  of  the  Delta  IV-H  rocket  vehicle.  The  Delta  IV-H,  as  designed,  is  dependent  on 
chemical  combustion  resulting  in  specific  impulse  values  for  its  1st  stage  engine  (the  RS- 
68)  of  420  s  and  for  its  2nd  stage  engine  (the  RL10B-2)  of  462.4  s.  The  goal  was  to 
identify  increases  in  performance  possible  with  alternative  power  sources  taking  into 
account  the  likely  increase  in  engine  mass  due  to  the  inclusion  of  shielding  and  changes 
in  materials.  Nuclear  fission  reactors  have  had  to  battle  environmentalists  and  safety 
concerns  over  the  course  of  their  existence,  leading  scientists  and  engineers  to  look  for 
propulsion  sources  capable  of  similar  perfonnance  without  the  threat  of  harmful 
radioactive  decay  products.  The  isomer  1 7*Hf”2  is  capable  of  high  levels  of  energy 
storage  (1.3  GJ/g),  releasing  only  gamma- radiation  during  its  decay.  The  absorption  of 
gamma-radiation  within  materials  can  lead  to  the  creation  of  charged  particles,  but  does 
not  lead  to  the  fonnation  of  the  long  lived  decay  products  common  to  the  fission  reaction. 
Based  on  the  amount  of  energy  stored,  the  length  of  time  for  which  it  is  stored,  and  the 
resulting  spectrum  of  decay,  this  material  is  viewed  to  have  the  potential  for  being  a 
source  of  energy  in  space.  This  work  examines  the  manner  in  which  this  material  could 
be  employed  within  a  rocket  to  facilitate  a  chain  reaction  within  the  material  and  achieve 
the  highest  levels  of  heat  transfer  to  a  propellant.  Due  to  its  similarity  to  nuclear  heat 
sources,  the  configurations  examined  were  the  NERVA,  PBR,  and  CERMET. 
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5.1  Conclusions 


1.  Application  of  a  nuclear  or  isomer  core  to  a  1st  stage  rocket  design  will  result 
in  a  significant  reduction  in  the  mass  savings  achieved  by  high  levels  of  specific  impulse. 
This  is  due  to  the  requirement  that  shielding  be  placed  around  the  entire  core  while  active 
within  the  lower  atmosphere.  An  isomer  design  without  shielding  will  be  dangerous  to 
humans  even  at  distances  of  miles  from  the  source,  and  there  is  potential  for  radioactive 
decay  products  to  be  expelled  in  the  rocket  exhaust  from  fission  core  designs.  The 
potential  hazards  outweigh  the  benefits  of  such  a  source  for  an  application  such  as  this 
since  no  fuel  savings  or  increases  in  payload  mass  are  predicted. 

2.  Based  on  reasonable  weights  for  shadow  shield  designs,  “light-weight”  core 
designs,  and  increased  levels  of  Isp,  achievable  with  hydrogen  flowing  through  the  fission 
designs,  significant  mass  savings  result  for  a  Delta  IV-H  2nd  stage  vehicle  (Tables  10). 
Systems  with  such  sources  can  capitalize  on  this  savings  by  carrying  more  payload  mass, 
achieving  greater  changes  in  velocity  (Appendix  F),  traveling  further  into  space,  or 
becoming  more  rugged  in  design. 

3.  A  PBR  dependent  on  the  fission  of  uranium- carbide  fuel  pellets  can  provide 
the  mass  savings  as  high  as  9,800  kg  when  applied  to  the  2nd  stage  mission.  Of  the 
fission  reactors  examined  in  this  study  the  PBR  has  the  characteristics  best  suited  to  take 
full  advantage  of  an  isomer  energy  source.  The  particle  bed  reactor’s  pellets  offer  large 
surface  areas  for  heat  transfer,  and  the  ability  to  manufacture  the  isomer  pellets  separate 
from  the  remainder  of  the  core.  Provided  the  pellets  can  be  swapped  out  of  the  core, 
potential  also  exists  for  a  reusable  engine. 
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4.  For  the  design  of  an  isomer  core  in  a  particle  bed  configuration,  high  levels  of 
efficiency  can  likely  be  obtained  if  core  masses  are  increased  to  specifically  absorb  the 
spectrum  of  radiation  being  given  off  as  a  result  of  the  decay.  Absorbing  gamma- 
radiation  requires  elements  high  in  atomic  mass  and  the  core  environment  requires 
materials  high  in  melting  temperature.  The  isomer  itself  will  have  a  melting  temperature 
of  2,506  K  in  the  metallic  fonn,  and  a  melting  temperature  of  4,173  K  if  hafnium  carbide 
is  used.  The  design  dependent  on  triggered  isomer  decay  for  its  energy  would  likely 
yield  somewhat  less  of  a  mass  savings  (in  the  range  of  6,000  to  9000  kg)  due  to  increased 
material  masses  and  lower  melting  temperatures  (in  the  case  of  metallic  hafnium),  but 
such  a  design  would  benefit  from  the  lighter  shield  mass  needed  to  protect  equipment 
from  its  spectrum  of  radiation  when  compared  against  a  fission  core.  Shield  loading 
factors  for  fission  designs  can  be  up  to  three  times  that  of  those  for  a  gamma  producing 
source  alone  (Appendix  H). 

5.2  Recommendations 

1.  While  the  triggering  of  1 7xHfTl2  is  still  an  area  of  active  research,  the 
investigation  of  photon  emission  from  various  materials  could  be  the  key  to  the 
successful  application  of  this  technology.  Ultimately,  achieving  a  chain  reaction  will 
depend  on  the  sensitivity  of  the  isomer  to  triggering,  but  equally  important  for 
applications  such  as  this,  are  means  for  creating  and  maintaining  the  population  of  trigger 
photons  within  the  core.  The  manner  in  which  this  is  carried  out  is  an  area  for  continuing 
research. 


85 


2.  The  investigation  of  tailoring  a  PBR  design  to  support  the  decay  of  the  isomer, 
in  addition  to  promoting  a  chain  reaction  within  the  included  isomer  material,  should  also 
include  an  analysis  of  the  material  quantities  needed  to  effectively  absorb  the  radiant 
energy  within  the  core  and  the  surface  area  requirements  for  convective  heating.  This 
will  have  a  large  impact  on  the  core  mass,  which  as  estimated  in  this  study,  assumes  a 
material  configuration  similar  to  that  designed  to  support  the  fission  reaction. 

3.  If  high  efficiencies  are  not  obtainable  as  predicted  by  the  estimates  in 
Appendix  I,  or  use  is  extended  beyond  unmanned  mission,  then  shielding  will  become  a 
more  important  aspect  of  this  design.  In  order  to  examine  weights  for  shield 
configurations  other  than  that  of  a  shadow  shield,  or  with  materials  other  than  lead,  build 
up  factors  specific  to  shield  shape,  type  of  source,  and  material  are  needed.  Only  a 
limited  number  of  potential  shield  materials  were  listed  in  the  sources  used  in  this  study, 
and  further  investigation  in  this  area  could  lead  to  the  optimization  of  shield  designs  for 
this  application. 

4.  Though  not  an  area  of  investigation  in  this  thesis  work,  an  argument  can  be 
made  for  all  of  these  systems  (nuclear  and  isomer)  as  sources  of  power  for  nuclear 
electric  systems  powering  deep  space  mission.  While  the  advantage  of  high  specific 
impulse  systems  is  more  clear  cut  in  the  deep  space  mission  category  due  to  chemical 
systems  inability  to  perfonn  such  missions,  this  does  not  mean  that  this  is  the  best  place 
for  a  system  such  as  those  considered  here  to  make  their  first  mark  on  the  U.S.  space 
program.  Shorter  missions  such  as  orbit  transfer,  and  possibly  a  trek  to  the  moon  would 
help  prove  the  reliability  of  these  systems  and  allow  engineers  to  perfect  their  designs 
before  traveling  into  the  outer  reaches  of  space  on  missions  spanning  decades.  An  isomer 
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core  could  even  be  returned  home  to  study  the  effects  of  its  operation  in  space,  since 
radioactive  decay  products  are  not  a  concern  with  this  type  of  energy  source. 

5.  The  isomer  1 7xHlm2  is  not  the  only  material  that  releases  energy  in  the  fonn  of 
high  energy  photons,  such  as  gamma-rays.  These  sources  of  energy  will  require 
materials,  such  as  those  considered  in  this  study,  to  capture  the  heat  of  the  radiation  if 
they  are  to  be  used  for  space  propulsion  applications.  Should  scientists  be  unable  to 
release  the  energy  stored  in  178Hfl2  reliably  and  on  demand,  then  work  such  as  this 
should  be  continued  specific  to  more  probable  fonns  of  energy  generation. 
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Appendix  A 


Material  properties  are  a  very  important  to  system  design  and  sizing.  Values  for 
density  (  p  ),  melting  temperature  (Tm),  and  thermal  conductivity  (k)  at  300  K  are  listed. 


Table  18.  Element  Properties 


Element 

Description 

Tm  (K)  (54) 

k  (W/cm-K)  at 
300K  (55) 

H 

Hydrogen 

0.08988 

- 

0.001815 

Air 

Main  Group 

0.001229 

- 

0.00024 

Be 

Alkaline  Earth 
Metal 

1.848 

1560 

2.01 

LiH2 

Alkali  & 
Hydrogen 

0.500 

962 

0.0005-0.003 

C 

Main  Group 

2.267 

3,800 

1.290 

A1 

Main  Group 

2.700 

933 

2.37 

B 

Metalloid 

2.3 

2,573 

0.274* 

Ge 

Metalloid 

22.5 

1,210 

0.64 

GaAr 

Semiconductor 

5.316* 

1,510* 

0.500 

Zr 

Transition  Metal 

6.511 

2,128 

0.227 

ZrC 

Ceramic 

6.7* 

3,813* 

0.210* 

Nb 

Transition  Metal 

8.570 

2,750 

0.537 

NbC 

Ceramic 

7.900* 

4,033* 

0.130* 

Mo 

Transition  Metal 

10.30 

2,896 

1.380 

Sm 

Lanthanide 

7.353 

1,345 

0.133 

Hf 

Transition  Metal 

13.31 

2,506 

0.230 

HfC 

Ceramic 

12.2 

4,173 

? 

Ta 

Transition  Metal 

16.65 

3,290 

0.575 

TaC 

Ceramic 

15.00* 

4,153* 

0.220* 

W 

Transition  Metal 

18.82 

3,683 

1.740 

wc 

Ceramic 

15.70* 

3,143* 

0.42 

Re 

Transition  Metal 

21.02 

3,443 

0.479 

Pt 

Transition  Metal 

21.09 

2,041 

0.716 

Os 

Transition  Metal 

22.5 

3,323 

0.9167 

Pb 

Main  Group 

11.37 

643 

0.353 

U 

Actinide 

18.74 

386 

0.276 

UC2 

Nuclear  Fuel 

11.3* 

2,623* 

3.3472 

UC 

Nuclear  Fuel 

13.5 

2.3849  (avg.) 

U02 

Nuclear  Fuel 

10.97 

2,827 

0.270 

*  denotes  items  found  on  Matweb  (56). 
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Appendix  B 


Table  19.  Mass  Attenuation  Coefficient  (|l/p  )  in  cnr/g 


Gamma-ray 

Energy 

(keV) 

10 

too 

200 

300 

400 

500 

600 

H 

0.3854 

0.2944 

0.2429 

0.2112 

0.1893 

0.1729 

0.1599 

Air 

5.120 

0.1541 

0.1233 

0.1067 

0.09549 

0.08712 

0.08055 

Be 

0.6466 

0.1328 

0.1089 

0.0946 

0.08471 

0.07739 

0.07155 

LiTC* 

0.3449 

0.1499 

0.1234 

0.1072 

0.09603 

0.08769 

0.08111 

C 

2.373 

0.1514 

0.1229 

0.1066 

0.09546 

0.08715 

0.08058 

A1 

26.23 

0.1704 

0.1223 

0.1042 

0.09276 

0.08445 

0.07802 

B 

1.255 

0.1391 

0.1136 

0.0986 

0.08834 

0.08065 

0.07460 

Ge 

37.42 

0.5550 

0.1661 

0.1131 

0.09337 

0.08212 

0.07452 

Ga 

34.21 

0.5197 

0.1619 

0.1123 

0.09325 

0.08236 

0.07487 

GaAr 

37.80 

0.5598 

0.1671 

0.1137 

0.09371 

0.08248 

0.07484 

Zr 

74.17 

0.9658 

0.2237 

0.1318 

0.10180 

0.08693 

0.07756 

Nb 

80.38 

1.0370 

0.2344 

0.1357 

0.10400 

0.08831 

0.07858 

Mo 

85.76 

1.0960 

0.2423 

0.1379 

0.10470 

0.08848 

0.07851 

Sm 

249.9 

2.9010 

0.5192 

0.2296 

0.14660 

0.11120 

0.09218 

Hf 

230.1 

4.1540 

0.7339 

0.3054 

0.18340 

0.13240 

0.10580 

Ta 

237.9 

4.3020 

0.7598 

0.3149 

0.18810 

0.13520 

0.10760 

W 

96.91 

4.4380 

0.7844 

0.3238 

0.19250 

0.13780 

0.10930 

Re 

101.1 

4.5870 

0.8119 

0.3339 

0.19760 

0.14090 

0.11400 

Os 

104.5 

4.6960 

0.8327 

0.3414 

0.20110 

0.14280 

0.11250 

Pt 

113.2 

4.9930 

0.8896 

0.3625 

0.21180 

0.14920 

0.11680 

Pb 

130.6 

5.5490 

0.9985 

0.4031 

0.23230 

0.16140 

0.12480 

U 

179.1 

1.9540 

1.2980 

0.5192 

0.29220 

0.19760 

0.14900 

*  weighted  average 

Mass  attenuation  and  mass  energy  absorption  coefficients  were  obtained  online  at  the 
National  Institute  of  Standards  and  Technology  (NIST)  website  (57).  The  gamma-ray 


89 


energies  chosen  where  based  on  the  spectrum  of  decay  emitted  from  1 78Hf”2  which  range 
from  12.7  to  574  keV  in  energy  (Fig.  1).  The  materials  listed  in  Tables  19  and  20  are 
current  materials  used  in  fission  reactors,  mediums  through  which  the  produced  radiation 
will  pass,  or  probably  materials  for  constructing  a  TIC. 

Table  20.  Mass  Absorption  Coefficient  (  (J-4/p  )incnf/g 


Gamma-ray 

Energy 

(keV) 

10 

too 

200 

300 

400 

500 

600 

H 

0.00985 

0.04063 

0.05254 

0.05695 

0.0586 

0.059 

0.05875 

Air 

4.742 

0.02325 

0.02672 

0.02872 

0.02949 

0.02966 

0.02953 

Be 

0.4225 

0.01838 

0.02353 

0.02548 

0.02620 

0.02639 

0.02627 

LiH2 

0.1222 

0.02066 

0.02666 

0.02889 

0.02971 

0.02991 

0.02980 

C 

2.078 

0.02147 

0.02655 

0.02870 

0.02950 

0.02969 

0.02956 

A1 

25.43 

0.03794 

0.02745 

0.02816 

0.02862 

0.02868 

0.02851 

B 

1.006 

0.0194 

0.02453 

0.02654 

0.02731 

0.02749 

0.02737 

Ge 

35.64 

0.3803 

0.06865 

0.03891 

0.03193 

0.02930 

0.02790 

Ga 

32.50 

0.34970 

0.06463 

0.03782 

0.03156 

0.02920 

0.02793 

GaAr 

36.01 

0.38380 

0.06921 

0.03916 

0.03210 

0.02943 

0.02802 

Zr 

71.50 

0.70800 

0.11640 

0.05420 

0.03885 

0.03311 

0.03025 

Nb 

77.54 

0.76080 

0.12470 

0.05705 

0.04026 

0.03396 

0.03085 

Mo 

82.75 

0.80420 

0.13160 

0.05919 

0.04117 

0.03437 

0.03104 

Sm 

227.8 

1.78600 

0.33660 

0.13260 

0.07620 

0.05411 

0.04334 

Hf 

197.5 

2.07500 

0.46450 

0.18530 

0.10350 

0.07044 

0.05409 

Ta 

202.1 

2.09200 

0.47840 

0.19150 

0.10690 

0.07248 

0.05545 

W 

92.04 

2.10000 

0.49130 

0.19730 

0.11000 

0.07440 

0.05673 

Re 

96.10 

2.10700 

0.50540 

0.20380 

0.11350 

0.07658 

0.05822 

Os 

99.40 

2.09200 

0.51500 

0.20850 

0.16100 

0.07813 

0.05923 

Pt 

107.8 

2.08100 

0.54130 

0.22160 

0.12330 

0.08265 

0.06230 

Pb 

124.7 

1.97600 

0.58700 

0.24550 

0.13700 

0.09128 

0.06819 

U 

171.1 

1.50200 

0.67460 

0.30500 

0.17320 

0.11520 

0.08494 
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The  expressions  1/p,  and  l/(l4  represent  the  average  distance  traveled,  measured  in 
centimeters,  before  a  photon  respectively  either  interacts  with  or  is  absorbed  in  the 
material  it  is  traversing.  The  shaded  blocks  indicate  the  materials  which  exhibit  the  best 
attenuation  or  energy  absorption  characteristic,  and  should  be  considered  first  as 
candidates  for  a  TIC  design.  In  order  to  get  some  idea  of  the  distance  a  photon  will  travel 
through  a  material  before  either  interacting  or  depositing  its  energy,  the  mean  free  path 
can  be  expressed  as  l/|l  or  1/p. A  measured  in  cm.  This  is  an  average  distance  value 
calculated  for  a  couple  of  different  energy  photons  as  shown  in  Table  21. 


Table  21.  Mean  Free  Path  for  Attenuation  and  Energy  Absorption 


Photon  Energy 

10-keV 

300-keV 

600-keV 

Element 

1/|I 

(cm) 

1 

(cm) 

1/|J, 

(cm) 

(cm) 

1/|J, 

(cm) 

(cm) 

Hf 

13.31 

0.000327 

0.000380 

0.246 

0.406 

0.710 

1.39 

Ta 

16.65 

0.000252 

0.000297 

0.191 

0.314 

0.558 

1.08 

Pb 

11.37 

0.000673 

0.000705 

0.218 

0.288 

0.705 

1.04 

W 

18.82 

0.000548 

0.000577 

0.164 

0.269 

0.486 

0.937 

Re 

21.02 

0.000471 

0.000495 

0.142 

0.233 

0.417 

0.817 

Pt 

21.09 

0.000419 

0.000440 

0.131 

0.214 

0.406 

0.761 

Os 

22.5 

0.000425 

0.000447 

0.130 

0.213 

0.395 

0.750 

h2 

0.08988 

28.9 

1,130 

52.63 

200 

6.944 

189 

Air 

0.001229 

159 

172 

7,626 

28,331 

10,101 

27,554 
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Appendix  C 


Specific  heat  equations  were  obtained  from  two  different  sources.  As  shown 
below,  the  equations  with  inputs  temperature  ( T)  and  the  molecular  mass  ( MM)  of  the  gas 
will  yield  the  specific  heat  value  ( cp )  of  the  indicated  gas  in  J/kg-K. 

From  Humble,  Henry,  and  Larson  (3:460): 


H?.  1000J/L/ 


'2'  ,,  =. 


MMHkg/  kmol 


[56.505  -  702. 74(—)“0-75  + 1165(— )_1  -560.7(—y15]kJ/kmol-K 


CH4:  c  = 


10007 /kJ 


100 

T 


100 

T 


100 
T 


- f— 67287-l-43974( — —  2487f( — 0  + 3238  ft — -)  ]kJ/  kmol-K 

MMCH  kg/  kmol  100  10 0  100 


C02:  c  — mOJIkJ — [_3i7357  +  30.529(— )05  -  4. 1034 (— )  +  0.024198 (— )2 ]kJ/ kmol- K 

p  MM co  kg! kmol  100  100  100 

From  Chemkin  data  file  (47).  Value  in  brackets  is  unites  and  Ru  is  in  J/kmol-K  here 

only. 

CO:  _  R'< 


MM,, 


-[3.025078  +  0.0014426897  -  0.05630828vl0^r  +  0.01018581xl0~s73  -  0.06910952x1  O^'T4 


N2:  =— —  [2.926640  +  0.00144879777-0.0568476M0~572+0.0109704xl0“s73-0.0675335M0~1374] 

p  mmN2 

°2'  q  —  — — [3.697578+  0.0006 135 1 977  —  0.01258842x10  5 72  +  0.0 177528 LxlO  y73  -0.01 136435x1 0“l3r4] 
p  MMn 


C,H«:  Ru 


MM, 


-[7.52521 7+ 0.01 8890347 -0.06283924xl0~472  + 0.09 1793  73x1 0~s73  -  0.048 1241 0x1 0“1274] 
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Appendix  D 


Particle-Bed  Reactor  sizing  equations  (3 1486)  were  obtained  for  three  reactor 
configurations  (7,  19,  and  37)  which  insure  criticality  in  a  fission  reactor  and 
geometrically  fit  together  to  form  a  nearly  circular  shape.  The  input  to  these  equations  is 
the  desired  core  power  ( Preq )  and  the  output  is  core  radius  ( RCOre )  and  core  height  (Hcore). 
For  7  Elements: 

Rcore  =  9.0958(1  Or'X  -1.3261(10)^  +  7.1 665(1 0)"4^ 
-0.1735/^+47.625 

Hcore  =  -0.000283/^  +  0.5203 Preq  +  26.06 

For  19  Elements: 

Rcore  =  -2.655(10 y'2Pr5eq  +  8.946(10 )~9  Pr4eq  - 1.1 703(1 0y5P?eq 
+  7.427(1 0)-3 Plq  -2.2955 Preq  +313.34 

Hcore  =  -4.027(1 0)-5/^  +  0.1427/^  + 17.9883 

For  37  Elements: 

RCOre  =  4.905(10) -2.881(10)-7P4  +  6.2522(1 0)"4^ 

-  0.5992/^  +  252.28 

Hcore  =  -6.502(1 0)-6/^  +  0.05009/^  +18.335 
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Appendix  E 


Figure  9  in  the  text  shows  that  radiation  dose  falls  off  according  to  the  inverse 
square  law  (1/r2),  and  the  data  in  Table  22  was  used  to  build  that  figure.  The  power  being 
produced  was  for  a  1st  stage  application  and  measured  41,794-MW.  The  cases  of  5%  and 
50%  of  the  produced  power  radiated  were  examined. 


Table  22.  Radiation  Dose  as  a  Function  o: 

f  Distance  from  the  Source 

Percent 

Radiated 

0.05  percent  radiated 

0.5  percent  radiated 

Distance  (m) 

Dose  rate  (rad/s) 

Log  scale 

Dose  rate 
(rad/s) 

Log  scale 

3 

4988800 

6.697996094 

49888000 

7.697996094 

10 

448990 

5.652236668 

4489900 

6.652236668 

20 

112250 

5.05018635 

1122500 

6.05018635 

50 

17960 

4.254306332 

179600 

5.254306332 

100 

4490 

3.652246341 

44899 

4.652236668 

200 

1123 

3.050379756 

11225 

4.05018635 

500 

180 

2.255272505 

1796 

3.254306332 

1000 

45 

1.653212514 

450 

2.653212514 

1609.344 

17 

1.230448921 

173 

2.238046103 

4828.032 

1.9 

0.278753601 

19 

1.278753601 

16093.44 

0.17 

-0.769551079 

1.7 

0.230448921 
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Appendix  F 


Calculating  the  velocity  change  requirements  for  the  two  stages  of  the  Delta  IV-H 
launch  vehicle  was  accomplished  in  the  following  manner. 

Stage  1: 

Orbital  Altitude  =  110  km 


v  = 


1 3. 986(10)5 km3  Is2 
V  6378km  +  1 1 0km 


1 .838km!  s 


Stage  2: 

Vehicle  mass  (m,)  =  41,533  kg 
Bum  time  (tb)  =  1,094  s  (For  GTO) 

Thrust  (F)  =  1 10,000  N 
Specific  impulse  ( Isp )  =  462.4  s 

Exhaust  velocity: 

ve  =  Ispg0  =  462.4s(9. 8lm/s2)  =  4, 536.144m/s 


Mass  flow  rate: 


110,000V 
4,536.144 ml  s 


24.2496 kg  Is 


Mass  of  propellant  used: 

m  p  =  m(tb)  =  24.2496kg/s(l,094s)  =  26,529.14kg 


Final  vehicle  mass: 

mf  =mi  —  mp  =  41,533kg  -  26,529.14kg  =  15,003.86kg 


Mission  velocity  change: 

mf  15  003  86kg 

Av  =  -v  ln( — — )  =  -4,536. 144m / s ln(  ’  ^)  =  4,6 1 8.6 lm / s 

m,.  41,533kg 


Maximum  velocity  change  possible: 

m,  14  333kg 

Av  =  -v  ln( — — )  = -4,536. 144m/ sln(  ’  2-)  =  4,826.1  bn/s 

m,  41,533kg 
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An  example  of  the  velocity  change  ( Av  )  effect  on  required  propellant  and  bum 
tune  is  shown  in  the  Tables  24,  25,  and  26.  Single  CERMET,  PBR  (19-Element),  or  TIC 
(19-Element)  used  without  shielding. 

Table  24.  Increases  in  2nd  Stage  Av  with  a  CERMET  Reactor  (No  Shielding) 


Av  (m/s) 

Hlextra  (kg) 

Burn  Time  (s) 

2,000 

18,138 

676 

3,000 

14,601 

961 

4,000 

11,442 

1,215 

4,619 

9,658 

1,359 

5,000 

8,620 

1,443 

6,000 

6,099 

1,646 

7,000 

3,847 

1,827 

8,000 

1,836 

1,989 

9,000 

39 

2,134 

Table  25.  Increases  in  2nd  Stage  Av  with  a  PBR  (No  Shielding) 


Av  (m/s) 

^extra  (^§) 

Burn  Time  (s) 

2,000 

19,113 

686 

3,000 

15,918 

981 

4,000 

13,021 

1,249 

4,619 

11,366 

1,401 

5,000 

10,396 

1,491 

6,000 

8,017 

1,711 

7,000 

5,861 

1,911 

8,000 

3,907 

2,091 

9,000 

2,136 

2,255 

10,300 

78 

2,445 
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Table  26.  Increases  in  2nd  Stage  Av  with  a  TIC  at  2,400-K  (No  Shielding) 


Av  (m/s) 

mextra  (kg) 

Burn  Time  (s) 

2000 

17,952 

674 

3000 

14,352 

957 

4000 

11,145 

1,209 

4619 

9,340 

1,351 

5000 

8,289 

1,434 

6000 

5,745 

1,634 

7000 

3,479 

1,812 

8000 

1,460 

1,970 

8800 

5 

2,085 
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Appendix  G 


In  sizing  the  PBR,  a  choice  between  3  configurations  known  to  satisfy  core 
criticality  requirements  and  fit  together  in  a  near  circular  configuration  will  result  in 
different  core  volumes  and  masses.  The  designer  can  choose  amongst  these 
configurations  based  on  the  space  available  for  the  reactor  or  based  on  the  limiting  the 
mass  of  the  rocket  as  much  as  possible.  The  core  will  the  smallest  radius  will  result  in 
the  smallest  shadow  shield  areas  and  subsequently  the  lightest  shadow  shield,  but  the  core 
with  the  smallest  radius  is  not  always  the  lightest.  The  values  that  were  set  in  this 
analysis  are  the  dimension  of  the  pellets  and  the  percentage  of  the  core  volume  that  is 
occupied  by  the  pellets  and  moderator  material.  Table  27  displays  these  values. 


Table  27.  PBR  Set  Values 


Pellet 

Outer 

Diameter 

(pm) 

Pellet 

Inner 

Diameter 

(|U/H) 

Total  Pellet 
Volume  (m3) 

Volume  of 

Inner 

Material 

(m3) 

Coating 
Mass  (kg) 

Percent 

Pellet 

Percent 

Moderator 

500 

300 

3.453(10)-' 

8% 

58% 

For  the  cases  of  18%  and  33%  of  the  power  generated  being  lost  to  radiation  the 
core  dimensions  are  calculated  and  shown  in  Table  28  below.  The  estimate  of  core 
weight  will  change  for  a  TIC  design  based  on  the  replacement  of  the  inner  pellet  material 
and  moderator  material  with  various  materials  that  are  beneficial  to  the  core’s  ability  to 
absorb  gamma- radiation  and  possibly  produce  trigger  photons  of  the  desired  energy.  The 
masses  calculated  here  assume  all  of  the  pellets  and  moderator  material  are  constructed  of 
the  metals  listed.  In  actuality,  the  pellets  and  moderator  could  be  a  mix  of  different 
materials,  and  the  moderator  dimensions  could  be  reduced  since  there  is  no  longer  a  need 
to  moderate  neutrons  within  the  core. 
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Table  28.  2nd  Stage  PBR  Reactor  Sizing  Analysis 


Configuration 

7  Element 

19  Element 

37  Element 

Percent  of 
Power  Radiated 

18% 

33% 

18% 

33% 

18% 

33% 

Core  Radius 
(m) 

0.3119 

0.3562 

0.3558 

0.3438 

0.8078 

0.5834 

2.187 

2.4196 

0.8076 

0.9215 

0.4239 

0.4742 

Core  Volume 
(m3) 

0.6691 

0.9111 

0.3212 

0.3421 

0.7228 

0.5071 

Core  Mass  (kg) 

1,071 

1,458 

514 

547 

1,157 

811 

Lead  Shadow- 
Shield  Mass 
(kg) 

421 

574 

548 

537 

2,823 

1,540 

Power  Density 
(GW/m3) 

0.769 

0.695 

1.602 

1.850 

0.712 

1.248 

Pellet  Volume 
Occupied  (m3) 

0.0535 

0.0729 

0.0257 

0.0274 

0.0578 

0.0406 

Coating  Mass 
(kg) 

282 

385 

136 

144 

305 

214 

Hf  Pellet  Mass 
(kg) 

154 

210 

74 

79 

166 

117 

W  Pellet  Mass 
(kg) 

218 

296 

104 

111 

235 

165 

Ta  Pellet  Mass 
(kg) 

193 

262 

92 

98 

208 

146 

Re  Pellet  Mass 

(kg) 

243 

331 

117 

124 

263 

184 

Os  Pellet  Mass 
(kg) 

260 

354 

125 

133 

281 

197 

Hf  Moderator 
Mass  (kg) 

5,165 

7,034 

2,480 

2,641 

5,580 

3,915 

W  Moderator 
Mass  (kg) 

7,304 

9,945 

3,506 

3,734 

7,890 

5,535 

Ta  Moderator 
Mass  (kg) 

6,461 

8,798 

3,102 

3,304 

6,980 

4,897 

Re  Moderator 
Mass  (kg) 

8,157 

11,108 

3,916 

4,171 

8,812 

6,182 

Os  Moderator 
Mass  (kg) 

8,732 

11,890 

4,192 

4,464 

9,433 

6,618 
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H.l  Thermal  Efficiency  of  Fission  Designs 

The  nuclear  fission  core  thennal  efficiency  values  displayed  in  Table  8  are  calculated  in 
the  fashion  shown  below. 


Thennal  Efficiency  Calculation  (Fission  PBR): 

Tabulated  values:  Power=  1,945  MW,  F=  333,617  N,  T=  3,200  K,  £  =  125:1 
Calculated  Power  =  1,598  MW 

Thermal  Efficiency  =  1,598/1,945  =  82%  Eff  =  ^core 

^tabulated 


Thermal  Efficiency  Calculation  (Fission  CERMET): 

Tabulated  values:  Power  =  2,000  MW,  F  =  445,267  N,  T=  2,507  K,  £  =  120:1 
Calculated  Power  =  1,798  MW 
Thermal  Efficiency  =  1,798/2,000  =  90% 

Thermal  Efficiency  Calculation  (Fission  NERVA): 

Tabulated  values:  Power  =  1,570  MW,  F  =  334,061  N,  7=2,361  K,  £  =  100:1 
Calculated  Power  =  1,299  MW 
Thennal  Efficiency  =  1,299/1,570  =  83% 


H.2  PBR  Fission  Design 

An  example  of  the  analysis  of  a  PBR  operating  on  the  principals  of  nuclear  fission  is 

canied  out  below. 

H.2.1  Rocket  Nozzle  Analysis 
Inputs 

Mission:  2nd  stage  of  Delta  IV-H  rocket 
Propellant:  H2  (molecular  mass  2.016  kg/kmol) 

Propellant  temperature  (7%):  3,200  K 
Chamber  pressure  (P02)'  6,176,000  Pa 
Nozzle  expansion  ratio  ( £  ):  285:1 
Thrust  (F):  1 10,000  N 
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Specific  heat  of  hydrogen  at  3,200  K: 


c  =  l000J,U  [56.505-702.74(2r2£)-°”  +l,165(2l^2£)1  _56<X7(2^V' 5]kJ/kmol-K  =18,642  (J/kg-K) 

PH‘  100  100 


2.016 kg!  kmol 


100 


Ratio  of  specific  heats: 


Y  = 


18,642 (J/kg-K) 


c„— ^  18,642 (J,kg-K)-  8.314.51  (Jlkmol-K) 


■  =  1.284 


MM  H  “  2.016(A'g/ kmol) 

Acoustic  velocity  at  stagnation  temperature: 


A  8,3\4A5J  kmol- K  , 

a0  =  y - A  =  1.284( - )3200A/  =411 6.6m  /  s' 

°  \  '  MM H  0  \  2.016/tg  / kmol 


Characteristic  velocity: 

an  411 6.6m/s 


c  = 


Y+l  1284+1 

1.284( - — - f^12 

y+1  1.284  +  1 


=  5468  Am  /  5 


Iteratively  solve  for  Mach  number  at  the  nozzle  exit: 

9  i  904  _ ,  1284+1 

0  =  M3(285)-[(-  —  )(1+  •  2  )]  2(1  284h2  ->  M3  =  6.7625 


1.284  +  1 


2M, 


Ratio  of  exit  pressure  to  chamber  pressure  (also  referred  to  as  stagnation  pressure) 


P 

_ e 


i+— — -m.2 


l=y 


,+(L^zll(6.7625)i 


-1.11  ( 1 0) 


Specific  impulse: 

5,468.4/;;/ 5  (1.284) 


'sP  =- 


9.81  m/s2 


_  .  1.284+1  1.284-1 

( - - - )( - - - )"i~284-I  {1  —  (1  ^11(10)  4)  I-2S4-  } 

1.284  -1  1.284  +1 


=  1036.65 


Exhaust  velocity: 


=  I spgo  =  1036.65(9.81w  /s2)=  10,169.046m  Is 


Mass  flow  rate: 

F 


110,000/V 


m  =  ■ 


/  g  1036.65(9.81//;  Is2) 


=  10.817/^/5 


Power  needed  to  heat  propellant  to  desired  temperature: 


p  =- 


-(32  0  0  025  -  300025) 


90,6/  /  L[0 +  56.505 (3200 -300)- 

z.Olokg/ kmol  0.25 

+ 1, 165(1 00)(ln(3200) - ln(300)) -  560-700°)  (32oo-a5  -  300"a5 )\U!  kmol- A"  =51 6MW 
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Core  radius: 


Rcm  =  -2.655(10)-'2(|i|)s  +  8.946(10)-’ (^|)4  -  1.1703(10)-’ (|i|)! 

+  7.427(10)-’(^-)2-2.2955(— )+  313.34  =  34.43cm 
0.82  0.82 


Core  height: 

Hcore  =  -4.027(1 0)-5  (^-) 2  +  0.1427(— )  + 17.9883  =  91.84cm 
0.82  0.82 


Core  volume: 

Vcore  =KR2coreHCore  =  7r(.3443m)2(-9184m)  =  0.3420m3 
Core  mass: 

mcore  =  PcoreVcore  =  (1,600kg/  m 3).  3420m 3  541kg 

Shadow  mass:  (Si  is  shield  loading  which  is  3,500  kg/m2  for  a  fission  shield) 
mshield  =  S,Ashield  =  S,KRle  =  3,500 kg/ m27t  (.3443 m)2  =  1,303.44kg 


H.2.2  Mission  Details 

Av  Requirement:  4,618.61  m/s 

Initial  Mass  (mi):  41,533  kg 

Propellant  Mass  ( mprop ):  27,200  kg 

Payload  Mass  (mi):  10,800  kg 

Pressure  Drop  through  the  Core  (Pdrop)'-  5% 

Final  vehicle  mass: 

-Av  -4,6 18. 61m/  s 

mf=m,(e  )  =  41, 533kg(em69M6m,s  )  =  26,372.06kg 

Mass  of  propellant  used  for  specified  mission: 

mp  =  m;  —mf  =  41,533kg- 26,372.06kg  =  15,160.94kg 

Fuel  savings  due  to  high  Isp : 

fnexlra=mprop-mp  =  26,529.14kg  -  15,160.94kg  =  1 1,368.20kg 

Current  system  structural  mass  without  chemical  engine : 
ms  =  3,490kg  -  30 1kg  =  3,1 89kg 
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Inert  mass  fraction  if  all  savings  is  used  for  structural  redesign  to  include 
,  Qnexlra+ms)  =  (11,368 .20kg+  3,189  kg) 

Jinert  mp+(mextra+ms)  15,160.94kg+ (1 l,368.20kg+3,189kg) 


fission  core 
0.4898 


Inert  mass  fraction  if  all  savings  is  used  for  additional  payload: 

_  (ms  +  mcore  +  mshiM)  _  (3,189frg+  547  kg+ 1,303 .44A-g) 

'  mp  +  (ms  +mcore+mshieU)  15,160.94kg+  (3,1 89  kg+5Alkg+ 1,303 ,44kg) 


Check  for  feasibility: 

Av  -4,61 8.61  ml  s 

finer/ SpS°  =  1  ~  0.4898e m6.6s(9.Slm/! s' !)  =  Q.6899  >  0 


H.3  Core  Breakup 

300  \im  diameter  sphere  ->  Vpi  =  1.4137(10)'1 1  m3 
500  \im  diameter  sphere  Vpo  =  6.545(10)" 11  m3 

Mass  of  fuel  in  fission  PBR: 

m  fuel  =  ^Mcore  =  ^  547kg  =  182kg 

Mass  of  moderator  in  fission  PBR 

moderator  =  J M core  =  ^547kg  =  365kg 

Mass  of  single  pellet  including  coating  and  fuel: 

mpelle,  =  p UCi  Vpi  +  p ZrC(Vpo  -  Vpi)  =  1 1,300kg/  m3  (1 .4137 (10)"11 )  +  6,730 kg!  m 3 

(6.545  (10)_1 1  - 1 .4137  (lO)'1 1 )  =  5 .0508(10)'7  kg 

Number  of  pellets  within  fission  PBR: 

,r  mf*  182kg 

Nn  =  — - —  = - ^ —  =  360,335,682 

P  mpellet  5. 0508(10)-7  kg 

Volume  occupied  by  pellets  in  fission  PBR: 

y  -  NrV?°  -  360,335,682(6.545(i0)-')  _  Q 
"  0.88  0.88 

Volume  occupied  by  moderator  in  fission  PBR: 
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v 

mod 


m 


mod 


Be 


— 365kg  =0.1975 m3 

1,848 kg/m3 


Percent  of  core  that  is  pellets  by  volume 


P  = 

1  bvp 


F 


0.0268m3 

0.3420m3 


=  0.0784  -  8% 


Percent  of  core  that  is  moderator  by  volume 


V 


mod 


bvm 


F 


0.1975m-' 

0.3420m- 


:  0.5775  ~  58% 


H.4  TIC  Mass  Changes 

2nd  Stage  Mission 

For  514.63  MW  (82%  efficiency) 

Core  sized  at  Vcore  =  0.3212  m3. 

Percent  of  pellet  volume  that  is  core  material  equals  21.6%  (300  pm  diameter) 
Percent  of  pellet  volume  that  is  coating  material  equals  78.4%  (500  pm  diameter) 

Total  mass  of  pellets  within  the  core  if  all  are  metallic  hafnium: 

mHf  =  PhAJO-2  I  6)Fw„  =  1 3,3  1 0kg/ nP  (0.08)(0.2 1 6)(0.32 1 2m3 )  =  73.87/cg 

Total  sum  of  all  pellets  coating  material: 

mzrc  =  Pzrc-^>yp(0-284)Fcorg  =  6,730  A:g/  m3(0.08)(0.784)(0.3212m3)  =  135.58&g 

Mass  resulting  from  replacing  moderator  with  metallic  hafnium: 

mHfmod  =  P  HfPbvmV core  =  1 3,3 1 0A:g  /  m 3  (0.58)(0.32 12m 3 )  =  2,479.6^g 


H.4.1  Shielding  the  19-Element  TIC 

Rcore  =  0.3558  m 

Power  radiated  from  the  core: 


P  = 

1  rad 


0.82 


P  = 

core 


422(10)61F 

0.82 


-422(10)61F  =  92.63(1 0)6W 


Mass  energy-absorption  coefficient  (  p A  /  p  =  0.0027  m2/kg)  is  an  average  value  for  air 
(40:369) 

Range  (r  =  3  m)  from  source  to  sensitive  components. 
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Allowable  dose  rate  (  D  =0.01  rad/s) 

Dose  rate  without  shielding: 

p  m  92fiV]0}6W 

Da  =  -^L.rjL  =  '  y  L  (0.0027 m  /kg)=  2,21 1.378 J/kg  =  221,137.84/W Is 

4 nr  p  4k  (3mY 

Relaxation  length: 

o 

Using  Fig.  7,  the  Buildup  Factor  ( B )  was  found  to  be  2.7  with  a  relaxation  length  of  17.9. 

Lead  shield  thickness: 

px  17.9  . 

x  = - = - =  0. 1208m 

\xPb  148.2 

Shield  loading: 

S,  =  pPbx  =  11, 400kg/ m3(.  1208m)  =  1, 377.12kg/  m2 
Shadow  shield  mass: 

mshield  =  AshieldS,  =  %RleS,  =7t  (0.3558m)2(l,377.12kg/m2)  =  547.69kg 


H.4.2  Maximum  Pellet  Temperature  for  TIC 

Thennal  conductivity  (k=  1.9  W/m-K) 

Maximum  static  temperature  of  propellant  flow  (73  =  2,400  K) 

Thickness  of  zirconium  carbide  coating  (x  =  100  pm) 

Outer  diameter  of  pellet  (d0  =  500  pm) 

Power  needed  to  heat  H2  to  2,400  K  (PCore=  422  MW) 

The  mass  of  hafnium  (m///)  pellets  in  the  core  is  detennined  from  the  core  size  in 
beginning  of  H.4. 


Surface  area  of  pellet: 

d 


A 


am^y 


471  ( 


500(10) 


-6  , 


m 


)2  =7.854(10) 


-7  . 


m 


Mass  of  single  hafnium  pellet  without  the  coating: 

mHfs  =  Vpip  Hf  =  1-4137(1 0)_11m3(13,310kg/m3)  =1.882(1 0)“7  kg 


Number  of  hafnium  pellets  in  the  core:  (very  sensitive  to 


m 


Hf 


m 


Hfs 


73.87  kg 
1.882(10)~7  kg 


=  392,610,592 


rounding) 
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Power  absorbed  and  available  per  pellet  for  propellant  heating  (if  evenly  distributed!): 


_Pcore  _  422(10)6r 

Np  392,610,592 


Pellet  interior  temperature: 

r  =T,+^  =  2,400K+  1-075^(100(10  )-»,„) 

'  '  kA  \.9W / m  -if  (7. 854(10)  “7  m) 


2£12K 


The  thermal  conductivity  is  varied  to  determine  the  affect  on  the  required  fuel  surface 
temperature  for  heating.  The  values  are  calculated  using  Eq.  [44]  and  graphed  in  Fig.  10. 


Table  29.  Generic  Case:  Variation  of  Temperature  with  Thermal  Conductivity. 


Thermal  Conductivity 
(W/m-K) 

Fuel  Surface 
Temperature  (K) 

0.2 

3084 

0.5 

2674 

1 

2537 

1.9 

2472 

3 

2446 

4 

2434 

100 

2401 

H.5  Sizing  a  NERVA  Fission  Core 

Typical  NERVA  power  density  (Pcien  =  1-57  GW/m3) 
Typical  NERVA  density  (  pNERVA  =  2,300  kg/m3) 

Required  power  for  2nd  Stage  mission  ( PCOre  =417  MW) 
Core  efficiency  =  83% 


Required  core  power: 


,  _  Pcore  _  4 1 7(1 0) 6  TP 
req  0.83  0.83 


502.4 MW 
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Radiated  power: 

Prad  Preq  -  Pcore  =  502.4Mtf  -417 MW  =  85.4 MW 


Core  volume: 


V„. 


502.4(10 )6W 
1.57(10  )9W/m3 


0.3200 m3 


Approximate  Hcore  =  2Rcore 
Core  Radius: 


R 


core 


1 


0.3200m3 
2 n 


0.3707 m 


Core  Height: 

Hcore  =  IRcore  =  2(0.3707/»)  =0.7413/7 1 


Core  Mass: 

mcore  =  P  NERVAKore  =  2.3  OOA'g  / /»’’ ( 0 . 3  200/»’ )  =  736 kg 

Shadow  mass:  (Si  is  shield  loading  which  is  3,500  kg/m2  for  a  fission  shield) 
mshield  =  S,Ashield  =  S,TiR2core  =  3,500kg/ m27t  (.3707m)2  =1, 511kg 


H.6  Sizing  a  CERMET  Fission  Core 

Typical  CERMET  density  (  pCERMET  =  8,500  kg/m3) 

Required  power  for  Stage  II  mission  ( PCOre=  435  MW) 
Core  efficiency  =  90% 


Required  core  power: 

p  _  Pcore  _  435(10)6r 
req  0.90  0.90 


483.3 MW 


Radiated  power: 

Prad  =Pm -Pcore  =  483.3MW-435MW  =48.3  MW 


Core  radius: 

Rcore  =0.0034Preq  +  20.79  =  0.0034(483.3)  +  20.79  =  22.43cm  =  0.2243m 
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Core  height: 

Hcore  =  0.0067  Preq  +41.418  =  0. 0067(483. 3)  +  41. 418  =  44.66cm  =0.4466 m 
Core  volume: 

Vcore  =7t  R]oreHcore  =jc  (0.2243m)2 (0.4466  m)  =  0.0706m3 
Core  mass: 

m(:ore  =  P CERMElKore  =  8,500^g  1 17?  (0.0706m3)  =  600.  Mg 

Shadow  mass:  (Si  is  shield  loading  which  is  3,500  kg/m2  for  a  fission  shield) 
mshield=S,Ashield  =S,TiR^ore  =  3, 500kg  / m2n  (0.2243m)2  =  553.2kg 
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Appendix  I 


Outlined  here  is  the  process  for  estimating  the  core  efficiency  of  a  gamma  producing 
source.  As  noted  in  the  paper,  a  spherical  layer  of  material  is  placed  around  a  point 
source  generating  only  gamma  photons  of  500  keV. 

Source  power  =  600  MW 

Radius  to  outer  surface  of  core  R  =  0.3m  =  30  cm 
Linear  attenuation  coefficient  for  lead  p  =  1.83  cm"1 

Linear  absorption  coefficient  for  lead  p4  =  1.035  cm"1 
Linear  attenuation  coefficient  for  tungsten  p  =2.66  cm'1 
Linear  absorption  coefficient  for  tungsten  p4  =  1.44  cm"1 


Source  strength: 

S  =  600(10)6  J/s(6.2415(10)18cF  /  J)( 


1 


500  (10)3  eV  /  photon 


)  =  7.5(10  )21  photons  /  s 


Uncollided  photon  flux: 


r°(30,x)  = 


_  7.5(1 0)21  photons  I s(e^x) 


=  6.63(10 )17  e  -'  photons  /  cm' 


4n R- 


471  (3  Ocnz)" 


For  the  case  of  px  =  4 

r°(30.2J86)=  75(10):'f,;,°Mm/"(e"4) 

47t(30cm)2 


1.215(10 )16  photons  I  cm 3  -s 


Rate  of  heat  deposition  at  materials  outer  surface  per  unit  volume: 

77(30,2.1 86)  =  1 .062(1 0)”6  erg  IMeV{  10~7  J  /  erg)paEoBa  (px)r°  (30,2. 1 86) 

=  1.062(1 0)-uJ/MeV(\.035cm-l)(0.5MeV)(2A9)\.2\5(l0y6photon/cm3 -s 
=  1462.36 W/cm3 

Volume  of  material: 

V  -  —  n[(R  +  xY  -  R3]-  —  7T [( 3 0 cm  +  2.186 cm)3  -  (30 cm)3  ]  =  26,568 cm3 
3  3 

Rate  of  heat  deposition: 

Hfinal(R,x)  =  H(R,x)V  =  1462.36JT  /  cm3  (26,568cm3)  =38851980ffi  =3S.9MW 


Mass  of  material: 

m  =  pPbV  =  11.34g/cm3(26, 568cm3)  =  301, 281g  =  301A:g 
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Core  efficiency: 


Eff=  1 


radiated 


=  1- 


produced 


38.9  MW 
600MW 


=  93% 


Values  for  the  two  materials  are  tabulated  as  follows. 


Table  30.  Spherical  Shell  Values  for  Lead 


|J,X 

X 

(cm) 

Ba 

r  °(R,x) 
(photons/cm3-s) 

H(R,x) 

(W/cm3) 

m 

H 

(MW) 

m 

(kg) 

1 

0.546 

1.51 

2.4396(10)1 7 

20,246 

6,288 

127.3 

71.3 

2 

1.093 

1.80 

8,879 

12,817 

113.8 

145 

■ 

2.186 

2.19 

1,462 

26,568 

38.9 

301 

7 

3.825 

2.61 

87 

49,009 

■a 

556 

10 

5.464 

3.01 

5 

73,735 

0.4 

836 

15 

8.197 

3.52 

0.04 

120,343 

0.005 

1,365 

Table  31.  Spherical  Shell  Values  for  Tungsten 


|lx 

X 

(cm) 

Ba 

r  °(R,x) 
(photons/cm3-s) 

H(R,x) 

(W/cm3) 

H 

(MW) 

m 

(kg) 

1 

0.376 

1.64 

30,585 

4,306 

131.7 

83 

2 

0.752 

2.16 

14,862 

8,720 

129.6 

168 

■ 

1.50 

3.26 

3,029 

17,827 

54 

344 

7 

2.63 

5.43 

251 

34,428 

8.64 

664 

10 

3.76 

8.37 

19 

48,077 

0.93 

927 

15 

5.64 

15.17 

0.235 

76,530 

0.018 

1,477 

Due  to  the  rough  nature  of  this  estimation  the  core  efficiencies  generated  for  the  nuclear 
fission  cores  were  used  in  the  analysis.  These  numbers  were  merely  generated  to  provide 
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a  level  of  confidence  that  high  efficiencies  with  mass  quantities  of  materials,  not 
prohibitive  to  a  TIC  design,  are  possible  for  heating  due  gamma-ray  emission.  The 
efficiencies  across  the  range  of  thicknesses  used  are  shown  in  Table  32. 


Table  32.  Spherical  Core  Efficiencies 


Material 

Lead 

Tungsten 

|xr 

x  (cm) 

Eff 

x  (cm) 

Eff 

1 

0.546 

78.8% 

0.376 

78% 

2 

1.093 

81% 

0.752 

78.4% 

4 

2.186 

93.5% 

1.50 

91% 

7 

3.825 

99.3% 

2.63 

98.6% 

10 

5.464 

99.9% 

3.76 

99.8% 

15 

8.197 

99.99% 

5.64 

99.99% 
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